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ABSTRACT 



Analysis and design of linear feedback control systems by application 
of the basic Mitrovic method is well established. However* little has been 
done in applying the basic Mitrovic method to nonlinear systems. The objec- 
tive of this work was twofold. First* Mitrovic 0 & method as applied to 
systems with a single gain-variable nonlinearity predicts limit cycles with 
results comparable to existing analysis methods. Secondly* an analog ccmput 
er study was made on a system with two gain-variable nonlinearities and the 
results interpreted by Mitrovic® s method. 

The analog computer results were obtained by a Bonner Scientific Corpor 
tion analog computer* Model 3100 s at the U, S. Naval Postgraduate School* 
Monterey* California, The writer wishes to express his appreciation for the 
assistance and encouragement given him in this endeavor by Professor George 
Thaler of the U, S, Naval Postgraduate School, 
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PART I 



PREDICTION OF LIMIT CYCLES BY MITROVIC'S METHOD FOR 
CONTROL SYSTEMS WITH ONE GAIN VARIABLE NONLINEARITY 
Introduction. 

Mitrovic's method is easily applied to nonlinear control systems if 
the nonlinearity can be described by an input-output characteristic in 
which the output is an odd single-valued function of the input. Such a 
nonlinearity can then be expressed as an instantaneous variable gain 
with zero phase shift. This property allows the nonlinear effect to be 
expressed as a variable coefficient of the system characteristic equation. 
For the nonlinearity in the error channel, such a system can only have limit 
cycles if the linear system is unstable since the nonlinearily contributes 
no phase shift. 

Mitrovic's method determines absolute stability for a linear system, 
which has fixed coefficients in the characteristic equation, by assuming 
two of the coefficients to be variable. The method derives equations 
for two coefficients with S and of the S-plane as parameters. 

With the JU) axis <3 = 0) as the S-plane mapping contour, the curve of 
one coefficient versus the other coefficient is plotted using different 
values of CO ^ . If this curve encircles a point defined as the working 
point or M-point in a clockwise direction for increasing Lt)p^ , the 
system is stable. The working point is determined by the fixed values 
of the two coefficients which were assumed variable. 

For a nonlinearity, as described above, one or more coefficients of 
the characteristic equation are variable. This manifests itself on the 
variable coefficient plane as a locus of possible working point locations. 
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Thus a limit cycle for the nonlinear system exists for the conditions 
at the intersection of the working point locus and the stability (Q = 0) 
curve. A typical Mitrovic plot is shown in Figure I. The frequency of 
the limit cycle is read from the stability curve at the intersection. 

T,he amplitude is determined approximately from the gain of the non- 
linearity obtained from the working point locus. 




Figure I 

Coefficient Plane of Mitrovic^ Method 
The amplitude determined in this manner will necessarily be less 
than the actual amplitude. ’ The gain of the nonlinearity is defined as 
the ratio of the instantaneous output to the instantaneous input and the 
actual amplitude could be determined accurately if during a cycle the 
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minimum /maximum instantaneous gain could be determined. However, the 
non-linear gain determined from the Mitrovic plot is a H psuedo - average** 
gain over a complete cycle, since this "psuedo - average 59 gain is deter- 
mined from an "average" M-point position over a cycle. 

The shape of the stability curve on a coefficient plane depends on 
the two coefficients chosen and the order of the characteristic equation. 

If the coefficients of the S° and are chosen variable,, then the stability 
curve for a third order system is a straight line on the Bo versus 0^ 
plane, for a fourth order system the stability curve is a parabola, and 
so on. If different coefficients are chosen variable, then the stability 
curve has a different shape on the corresponding coefficient plane. 

The working point locus depends also on the linear system and speci- 
fically the number of zeros in the loop transfer function. For no zeros 
in the loop transfer function and the nonlinearity in the forward path, 
the working point locus is as straight line segment parallel to the 0^ 
axis since the variable nonlinear gain appears only in the constant co- 
efficient of the characteristic equation. For a linear system with one 
zero in the loop transfer function, the working point locus is again a 
straight line segment with a specific slope on the 0^ versus 0j plane. 
Other systems can be conceived with different placements of the nonlinear- 
ity but the basic working point locus is always a straight line segment 
on the appropriate coefficients plane. 

The use of Mitrovic* s method for a nonlinear system for which the 
linear loop transfer function has two or more zeros becomes more compli- 
cated. In this case the nonlinear gain appears in three or more co- 
efficients of the characteristic equation. With only two assumed variable 
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coef f icient s s the stability curve becomes a family of curves with the 
nonlinear gain a parameter. To determine the limit cycles for such a 
system the intersection of the straight line working point locus and 
the family of stability curves must be found. 

The various systems chosen to analyze in the chapters of Part One 
are completely arbitrary. The author had no specific system cr specific 
device in mind when the block diagram was drawn. Consequent ly 9 all signals 
in the various systems are referred to in volts. The following guidelines 
were followed in order to simplify the analysis without disrupting the 
arbitrary nature of the system* 

1. The poles and zeros of the linear loop transfer function 
were in general chosen one octave apart. This was done in order to 
sketch the actual Bode magnitude curve in the describing function anal- 
ysis from the asymptotic plot by applying the usual corrections at the 
corner frequency and two octaves either side of the corner frequency. 

2. The magnitudes of the poles, zeros, and gain constants 
were chosen to eliminate magnitude and time scaling procedures in the 
analog computer simulations of the system. 
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CHAPTER I 



SATURATION IN ERROR CHANNEL 



I-A Loop Transfer Function with No Zeros. 

i The block diagram of the system chosen to analyze is given in Figure 

’ 4 




’ Figure I-A 

Saturation in Error Channel 
Loop Transfer Function with No Zeros 

I-A- 1 . Mitrovic's Method. 

Assigning the variable gain N to the nonlinearity, the character- 
istic equation becomes, 

S 3 +BS 2 + ZS t Idhl =0 , (i-A-i.i) 

Assuming the last two coefficients variable the characteristic 
equation is rewritten, 

■S 3 '* - 3 5 + 0 ( S + 8a =0 (I-A-1.2) 
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From Appendix A, Mitrovic's equations for and are, 

8 0 = -[3t On<t>^) t O>n 3 0 2 (<S)] (i-A-1.3) 

and 

B,= 3u> n 0 2 (£) + U) n Z tf 3 (%) . (I-A-1.4) 

Substitution of the values of the 0 functions Appendix A for ^ = 0 
gives the parametric equations of the stability curve as 

8* =3 Wo 2 - (I-A-1.5) 



and 




(X-A-1.6) 



From comparison of equations I-A-l.l and I-A-1.2, it is seen that 
the equations describing the M-point locus are 



e 0 = i on 



(I-A-1.7) 



and 



B, 




(I-A-1.8) 



Thus the set of four parametric equations I-A-1.5 through I-A-1.8 
describe the two curves whose intersection is desired for limit cycle 
determinat ion. 

At the intersection, Equation I-A-1.5 equals Equation I-A-1.7 and 
Equation I-A-1.6 equals Equation I-A-1.8. Thus, 

I0 n Z = 2. (i-a-i.9) 
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and 



IQ N ave = 3i X> n z . (i-A-i.io) 

Solution of these two simultaneous equations yields, 

% 

UC^= 1.414 rcj/sec . , (i-a-i.h) 

and 

i 

Nave ~ 0*6. (i-A-i. i2) 

Since the nonlinear gain N is defined as the ratio of the output 
to the input, 

^ _ r 

£ flqa/ — \ | (I-A-1.13) 

Thus the frequency of the limit cycle is given by equation I-A-l.ll 
and the approximate amplitude by I-A-1.13 or 

= 8 . 33 Volts . (i-A-i.i4) 

The above solution, is presented graphically in Figure I-A-l.l. 




Graphical Solution for Saturation in Error Channel, Loop Transfer 

Function with no Zeros 
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The limit cycle is a stable limit cycle as determined by the follow- 
ing reasoning. If the input signal, 0^ , is large enough such that ft 

is small enough to place the instantaneous M-point at say tA ^ in Figure 
I-A-l.l. The system is instantaneously stable and the oscillations de- 
crease. As the oscillations decrease, N increases moving the instantan- 
eous M-point above the stability curve. Here the system is unstable and 
the oscillations increase thereby decreasing N and moving the M-point back 
into the stable region. Eventually a dynamic equilibrium will occur in 
which the "pseudo average" M-point lies on the intersection of the two 
curves. If the input signal is small enough, N will be large enough 
such that the instantaneous M-point lies at say at which point the 

system is unstable. By similar reasoning as above, the "pseudo - average" 
M-point will again eventually end up at the intersection. 
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I-A-2, Root Locus Method. 



The root locus method of analyzing a system with such a non-linearity 
consists of construction of the linear root locus and determining the 
point where the root locus crosses the J60 axis. 

The system characteristic equation is 

5 31 + 3s* + 2 s + ION = 0 . (i-a- 2. i) 

Substitution of $=j£e> and simplifying 

-Ju> 3 - 3 LO 2 +12UJ + ION =0 . U-A-2.2) 

By requiring that both the real and imaginary part of equation I -A- 
2.2 go to zero independently , 

-j U) 3 4 J 2 00 =0 (I-A-2.3) 

and 

- 3 co 2 + |0N - 0 , (i-a-2.4) 

Dividing equation I-A-2.3 by the frequency is obtained as 

U) =s 1.414 rad/sec. 

By substituting to into equation I-A-2.4 the ,9 pseudo - average 0 
gain over a complete cycle is obtained as 

Nai/e = 

from which the approximate amplitude of the limit cycle is determined 
as 

- 8 - 33 voits . 



9 



A sketch of the root locus solution is given in Figure* I-A-2. 1 . 




Root Locus Solution for Saturation in Error Channel 
Loop Transfer Function with No Zeros 
The root locus method also gives a stable limit cycle by the follow- 
ing reasoning. If the input signal is large enough, N will be small 
enough such that the instantaneous root location is at in Figure 

I-A-2.1. The system is stable and the oscillations decrease. As the 
oscillations decrease N increases moving the instantaneous root into the 
right half plane. The system is now unstable, the oscillations increase 
causing N to decrease, and the instantaneous root rrioves back into left 

t 

' ? * 
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half plane. Eventually a dynamic equilibrium is established in which the 
"average’’ root location is at the imaginary axis cross-over. If the 
input signal is small f N is large placing the instantaneous root in the 
right half plane at • By a similar reasoning the "average" root 

location over a cycle will again lie at the imaginary axis cross-over. 
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I-A-3. Describing Function Method. 

Application of describing function theory to the given system entails 
plotting the linear loop transfer function and the negative reciprocal of 
the describing function on the gain-phase plane (Nichols plot). Limit 
cycles of the system are then given by the intersections of the two curves 
on the gain-phase plane. 

The describing function for a saturating element with unity gain in 
the linear region is 

Gq - R -f R V I - R 2, J L Q . (i-a- 3 . 1 ) 



where 



- - 



R - — * (i-A-3.2: 

t max . 

The linear loop transfer function in frequency response form is 



6W>= J^OuniW^I) ' a ’ a - 3 - 3 

Since has no phase shift thenl/G^ has an associated phase angle 
of -180°, and the 1/G^ curve lies entirely on the -130° axis of the gain- 
phase plane. Consequent ly s the only region of the linear loop transfer 
function which is of interest is the region in which the phase shift is 
-180°. The values of gain in db and phase of the linear loop transfer 
function in the region of -180° which were used to plot on the gain- 
phase plane are given in Table I-A-3.1. These values were obtained from 
a Bode diagram plot of the linear loop transfer function. 
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TABLE I -A- 3.1 



to 


| G0^)ldb 


/6QU) d&3 


1.0 


9.8 


-162 


1.1 


8.5 


-167.5 


1.2 


7.1 


-172 


1.3 


5.8 


-176 


1.4 


4.5 


-179.5 


1.5 


3.3 


-183 


1.6 


2.2 


-186 



The values of 1/G^ in db were calculated for various values of R 
according to equation I-A-3.1 and appear in Table I-A-3.2. The gain- 
phase plane of the system is shown in Figure I-A-3.1. 



TABLE I-A-3.2 



R 


Go 


V So 


VSd ^ 


0.0 


0 


CO 


oC 


0.1 


0. 127 


7.88 


17.92 


0.2 


0.253 


3.95 


11.94 


0.3 


0.376 


2.66 


8.5 


0.4 


0.496 


2.02 


6.1 


0.5 


0.609 


1.642 


4.32 


0.6 


0.716 


1.396 


2.9 


0.49 


0.598 


1.673 


4.48 
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By interpolation of Figure I-A-3.1, the intersection occurs at 



(JO = 1.41 rad/sec 

and 



R = 0.495 # 

Thus from Equation I-A-3.2 the amplitude of the limit cycle is 



6 (BOX' 10 -* volts • 

Assuming the 1/G^ curve to be the locus of the critical point of 

the system, the describing function method predicts a stable limit cycle 

by the following reasoning. If the input signal* ©X is large R will 

be small and the instantaneous critical point lies above the curve 

say at point R in I-A-3.1. The system is momentarily stable due to a 
a 

positive phase margin and the oscillations will decrease. As the ascii- 
lations decrease, R increases and the critical point moves below the 
GOu)) curve. The system is new unstable due to a negative phase 
margin and the oscillations increase moving the critical point above the 
curve. Eventually a dynamic equilibrium is reached in which the 
"average" critical point over a cycle lies at the intersection of the two 
curves. If the input signal is small, R is large and the critical point 
is momentarily below the G£)U))curve at . The system is unstable 

and the oscillations increase. Similar reasoning shows again that the 
"average" critical point over a cycle will eventually lie at the inter* 
section of the two curves. 
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I-A-4. Analog Computer Simulation, 



The simulation of the given nonlinear system is given in Figure 
X-A-4.1. The coefficient pot settings with associated resistances and 
capacitances for real time and real magnitude scaling are given In Table 

I-A-4. 1. 



TABLE I-A-4.1 



Pot Setting 


Associated Elements 


E i 












a i = 100 


For 


Input Step of E, 


. volts i 

l 


ifid 




R i 


« 1 Meg, 


R fl 


M 1 


Meg 


a 2 = 1.0 


R 2 


= 1 Meg 


R fl 


« 1 


Meg 


a 3 = 1.0 


R 3 


= 1 Meg 


R f2 


« 1 


Meg 


a. = 1.0 
4 


R 4 


= 1 Meg 


C fl 


® 1 


uf 


a c - 1.0 

J 


R s 


~ 1 Meg 


C fl 


» 1 


uf 


a, = 1.0 
6 


R 6 


= 0.1 Meg 


C f2 


= 1 


uf 


a ? = 1.0 


R 7 


- 0.5 Meg 


C f2 


« 1 


uf 


a g = 1.0 


R 8 


= 1 Meg 


C f3 


= 1 


uf 


a = 0.05 


For 


a 5 Volt Saturation 






Initial 


Voltage 








Setting 













To accurately set the back-biasing pots, * 9 a" 9 a sine wave generator 
was connected to the input of the saturation simulator. The two pots 
were then adjusted until the saturated output voltage was exactly 5 volts 
on a recording device. The frequency of the generator was set at O.'Ol cps 
and the static characteristic of the saturation simulation was obtained 
by an X-Y recorder. This characteristic appears as Figure X-A-4.2. 

The system was first excited with a 2 volt step input signal. The 
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system was found to go in to the limit cycle shown in recorder trace 
I-A-4.1. The system was then excited with a 20 volt step input and was 
found to go into the same limit cycle {recorder trace T-A-4.2). From 
the recorder traces the limit cycle is a stable limit cycle with an 
angular frequency of 1.395 rad/sec and 10.0 volts amplitude. The phase 
portrait for the system appears in Figure I-A-4.3. 
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CHART NO. RA 2921 30 BUUSH TNSTKUMEXTS 0 \ 




Recorder Trace X-A-4,3 % 

Limit Cycle Measurement y Error Channel Saturation, Loop 
Transfer Function with $ o Zeros, • 



Scales: - 1 Volt/Line X ** 0,5 Volts/ Line 

Time - 0,2 Second/Division 
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I-A-5. Comparison of the Four Analytic Methods in Limit Cycle Pre- 
diction for the System of Figure X-A. 



Method 


Frequency 

(rad/sec) 


Amplitude 

(volts) 


Mitrovic 1 s 


1.414 


8.33 


Root Locus 


1.414 


8.33 


Describing Function 


1.41 


10.1 


Analog Simulation 


1.395 


10.0 
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I-B. Loop Transfer Function with One Zero. 

The block diagram of the system chosen to analyze is shown in 
Figure I-B. 




, Figure I-B 



Saturation in Error Channel Loop 
Transfer Function with One Zero 
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I-B-l. Mitrovic's Method. 



If N is the instantaneous gain of the nonlinearity the character- 
istic equation for the system is 

S 4 t I3S 3 + 44s 2 +(3Z + 300N)s + C>OON-0 , (i-b-i.d 

Assuming the last two coefficients to be variable, the characteris- 
tic equation becomes, 

S^ + i35 3 4*44s 1 + 8, s ■+ Bo = 0 (i-b- 1 . 2 ) 

where 

B, = 32 + 300 N (I-B-l. 3) 

and 

B 0 = 600 N , (I-B-l. A) 

From Appendix A, Mitrovic f s equations for and are, 

8 0 = -[44 + w„ 4 few)] (i-B-i. 5) 

and 

8, = 44Wn^(<0 + i3to + Loffyiy) . (i-B-i. 6> 

Substitution of the values of the 0 functions from Appendix A for 
= 0 gives, 

8 0 - 44 U> n 2 - (I-B-l. 7) 

and 

6, - I3(0 0 a . (i-B-i. 8) 

Thus the sets of parametric equations I-B-l. 3, I-B-l. 4, I-B-l. 7* 
and I-B-l. 8 describe the stability curve and the working point locus on 
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the B versus B plane. To find the intersection of two curves solve 
o 1 

Equation I-B-1.8 for ^ and substitute in Equation I-B-1.7 to obtain 



a - AA o __ JV 

- 1 3 



(I-B-l.S) 



Solving Equation I-B-1.4 for N and substituting in Equation I-B-1.3 
gives, 

8, - 32 + G.S’ B c (i-B-1.10) 

At the intersect ionj equations I-B-1.9 and I-B-1.10 must be satis- 
fied simultaneously. Solving Equation X-B-l. 10 for and substitut- 
ing in Equation I-B-1.9 yields, 

8, z -234 8 , - 10816 = o (i-B-i.il) 

Solving for gives 

B, = 273. 54 s and B. = -39.54 
1 1 

The negative value is meaningless, thus substitution of the posi- 
tive value into equations I-B-1.8 and X-B-1.3 gives the frequency and the 
"average" gain respectively, as 

CD?) “ 4.587 rad/sec 

and 



N av , e = °- 8051 

The approximate amplitude of the limit cycle is given by 

E 



max 



sat 

N 



or 



max 



5.0 

0.8051 



6.21 Volts 



The graphical representation of the solution is shown in Figure 

I-B-l.l. 
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Figure I-B-l.l 

Graphical Solution for Saturation in Error Channel, 
Loop Transfer Function with One Zero 



Mitrovic's method predicts a stable limit cycle by feli.ovi-.tg simii 



reasoning to that given in SectLcn I-A-l of this chapter. 
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I-B-2. Root Locus Method. 



The system characteristic Equation is 

s 4 + 13s 3 + 44s 2 + (32 + 300 N)s + 600 N = 0. 

Substitution of s = J LD and simplifying s yields 

u; 4 -Jl3LO a -4^liO i +^(32-f30DW)Lo+ bOOH - O 



(I-B-2.1) 



(I-B-2. 2) 



By requiring that both the real and imaginary parts of equation 
I-B-2.1 go to zero simultaneously s then 

600 N = 0 (I-B-2 . 3) 

and ~ 0)V 32U3-I- J 30060 N = 0. (I-B-2. 4) 

Dividing Equation I-B-2. 4 by JOJ and solving for N gives 

13 ,..' 2 32 



N 300^ 300 



# (I-B-2. 5) 

Substituting this into Equation I-B-2, 3 gives the following equa- 
tion for ID 

CO* -ISu/- 64 = 0 « (I-B-2. 6) 

2 

Solution of Equation I-B-2, 6 for 0J gives 

6U 2 = 21.042, 

from which 

CO = A, 587 rad/sec, 

2 . 

Substitution of (x) into Equation I-B-2. 5 gives the "average" non- 
linear gain over a cycle as 

N = 0.8051 
ave 

From this the approximate amplitude is obtained as 

£. , = n Wi = 6 - 21 Volts 

max 0.8051 

A pictorial representation of the root locus solution appears in 



Figure I-B-2.1. 
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The root locus method also predicts a stable limit cycle by follow 
ing similar reasoning to that given in Section I-A-l of this Chapter. 




Figure I-b-2.1 



Root Locus Solution for Saturation in Error Channel 
Loop Transfer Function with One Zero 

I 
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I-B-3. Describing Function Method. 

The describing function is the same as that given in Equation X-A-3.1 
of this Chapter. 

The linear loop transfer function in frequency response form is s 



a (M) = te.7f(0.5--)U>tl) 



(I-B-3.1) 



Jo)W+i)(o,is:Ju;+ixo,i2rjLu+/) * 

Table I-B-3.1 gives the magnitude and phase of Equation X-B-3,1 
which were used to plot the linear loop transfer function on the gain- 
phase plane of Figure I-B-3.1. These values were obtained from a Bode 
diagram plot of the linear loop transfer function. 



TABLE I-B-3.1 



0L> 


IGoortJjb 


/60«) deg 


4.05 


4 


-174.5 


4. 3 


3 


-177 


4.55 


2 


-180 


4.85 


1 


-183 


5.15 


0 


-190.5 



Table I-B-3.2 gives the values of 1/G^ in the area of interest. 
These values were obtained from Equation I-A-3.1. 

TABLE I-B-3.2 



R 


Gd 


1/ Go 


\/&D 


0.6 


0.716 


1.396 


2.9 


0.68 


0.794 


1.26 


2.01 


0.69 


0.803 


1.246 


1.92 


0.7 


0.812 


1.232 


1.83 
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By interpolation of Figure I-B-3.1, the intersection and thus the 
limit cycle occurs at 

L o =4.55 rad/sec. 

and 



R = 0.682 

From Equation I-A-3.2, the amplitude of the limit cycle is then 

£ „ = — = 7.33 Volts 

max 0,682 

The describing function predicts a stable limit cycle by following 
similar reasoning to that given in Section I-A-3 of this Chapter. 
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I-B-4. Analog Computer Simulation. 

The analog simulation of the given nonlinear system is shewn in 
Figure I-B-4.1. The coefficient pot settings with associated resistances 
and capacitances for real time and real magnitude scaling are given in 

Table I-B-4.1. 



TABLE I-B-4.1 



Pot Setting 


Associated Elements 


E i 




a i = Too 


For Input Step of Volts 




and = 1 Meg, = 1 Meg 


&2 ~ 1.0 


R 2 “ 1 Meg, R fl - 1 Meg 


a 3 = 1.0 


R 3 = 1 Meg s R f2 *= 1 Meg 


3 4 = l '° 


R^ - 0.5 Meg, s 1 uf 


= 0.4 


Rj = 0.1 Meg, = 1 uf 


a, = 1.0 
6 


R fi = 0.1 Meg, R < , 3 = 1 Meg 


a ? = 1.0 


R 7 = 0.1 Meg, R f3 = 1 Meg 


a 8 ' °' 3 


R q - 0.1 Keg, C,« - 1 uf 


« 9 - 1.0 


Rg = 1 Meg, C f2 = 1 uf 


a io " uo 


R 10 = 0.1 Meg, C f3 * 1 uf 


a H = 


R n = 0.1 Meg, C f3 = 1 uf 


a 12 - l '° 


R 12 - 1 Meg, C f4 - 1 uf 


a n . 1.0 


R^ 3 = 1 Meg, = 1 Meg 


a » 0.05 


For a 5 volt 


initial Setting 


Saturation Voltage 



To accurately set the back-bias pots "a" a sine wave generator was 
connected to the input of the saturation simulator. The two pots were 
then adjusted until the saturated output voltage was 5 volts on a re- 
cording device. 
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The static characteristic for the saturation simulator appears as 
Figure I-A-4,2 in Section I-A-4 of this Chapter. 

The system was first excited with a 4 volt . input signal* The 
system was found to go into the limit cycle shown in recorder trace 
I-B-4.1, The system was then excited with a 10 volt step input signal. 
The system was found to go into the same limit cycle (recorder trace 
I-B-4.2). From these two recorder traces, the limit cycle is seen to be 
stable with an angular frequency of 4.58 rad/sec and 7* 2 volts amplitude* 
The phase portrait for the system appears in Figure I-B-4.2, 
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I-B-5. Comparison of the Four Analytic Methods In Limit Cycle Prediction 



for the System Given in Figure I-B. 



Method 


Frequency 

rad/sec 


Amplitude 

(volts) 


Mitrovic 


4.587 


6.21 


Root Locus 


4.587 


6.21 


Describing Function 


4.55 


7.33 


Analog Simulation 


4.58 


7.2 
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CHAPTER II 



DEAD ZONE IN ERROR CHANNEL 
II- A Loop Transfer Function with No Zeros. 

/ 

The block diagram of the arbitrarily chosen system which is to be 

• 4 ' 

analyzed is given in 'Figure II-A. 




r- , t \ - ^>0 

W S(S+IXS+2XS+8) 



Figure Ii-A 

Dead Zone in Error Channel Loop 
Transfer Function with No Zeros 



II-A-1. Mitrovic's Method. 



The characteristic Equation of the given system is, 

s 4 +lls 3 + 26s% Its-f (oO l\l = o . 



(II-A-1. 1) 



Assuming the last two coefficients variable the characteristic equa- 
tion becomes, 

S 4 + II s 3 -h 2^S Z -f 8, S 4 6o = o (II-A-1. 2) 

where the parametric equations of the working point locus are* 

B = 60 N (II-A-1, 3) 

o 

and 



B 1 = 16 „ 

From Appendix A, Mitrovic's equations for and are s 



(II-A-1. 4) 



6 0 = (#)+ (II-A-1. 5) 



and 



6, = 2 L>u> n <h0£)+ l luffcL® + u n 3 fi 4 (2) • (ii-A-i. 6) 



Substitution of the values of the 0 functions from Appendix A for 
=0 gives the parametric equations of the stability curve as 

3 0 - 26(ja n 2 ~k>n 4 (h-a-1.7) 

and 

6,-= Ilu0 n 2 ‘ , (II-A-1. 8) 

At the intersection of the M-point locus and the stability curve, 
equations II-A-1, 3 and II-A-1. 4 and Equations II-A-1, 7 and II-A-1,8 
must be satisfied simultaneously. Thus at the intersection, 

n&yf" 16 (ii-a-i. 9) 



42 



from which the angular frequency of the limit cycle is determined to be, 

LO^ = 1.206 rad/ sec. 

Also at the intersection 

= 40N . (II-A-1.10) 

' Substitution of, UOt^ from Equation II-A-1.9 i.nto Equation II-A-1.10 
gives the "average" gain of the nonlinearity over the cycle. as, / 

N = 0.595 , 

ave 

To obtain the approximate amplitude of the limit cycle, reference is 
made to Figure II-A-1.1 which is a sketch of the nonlinearity static 
characteristic. 




Figure II-A-1.1 

Limit Cycle Amplitude Calculation from the 
Dead Zone Characteristic Knowing N flve 
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From the figure the following equations are obtained* 



and 



£ = 5 + c 

max 

c = a cos 45° = 0.707a 

5 Sin A 
a ~ Sin B 

A = Tan” 1 N 

ave 



(II-A-l.il) 

(II-A-1.12) 

(II-A-1.13) 

(II-A-1.14) 



B = 180 - A - (180 - 45) , (II-A-1.15) 

Entering equation II-A-1.14 with the "average" gain over a cycle* 

A = Tan” 1 0.595 = 30.8° . 

Substituting this in Equation II-A-1.15 gives 

B = 45° - 30.8° = 14. 2 Q # 

Substituting for A and B in Equation II-A-1.13 gives* 



a = 5 



sin 30.8 



= 10.469 



sin 14.2 

Substituting for "a" in Equation II-A-1.12 gives 

c = 0.707 (10.469) - 7.402. 



And finally Equation II-A-1.11 gives 

£ = 5 + 7.402 = 12.402 Volts 

max 

The graphical solution of the intersection is sketched in Figure 



II-A-1.2. 



Mitrovic's method predicts an unstable limit cycle by the following 
reasoning. If the input signal is small* N will be small and the in- 
stantaneous M-point will be in the stable region such as M b in Figure 
II-A-1.2. The system is stable and the oscillations will decrease in ampli- 
tude. As the oscillations decrease* N becomes even smaller driving the 
M-point even further into the stable region. Consequently the M-point 
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Figure II-A-1.2 

Graphical Solution for Dead Zone in Error Channel,- 
Loop Transfer Function with No Zeros 

can never go into the unstable region and all oscillations will eventual- 
ly die out. If the input signal is large, N is relatively large and the 
instantaneous M-point lies in the unstable region such as M . The system 

Si 

is unstable and the oscillations increase, N increases Coward N = 1.0, 
and the M-point moves farther into the unstable region. Consequently, 
the oscillations increase indefinitely in magnitude. 
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XI-A-2. Root Locus Method, 

Substitution of s = -> LO Into the characteristic equation given in 
Equation IX-A-1.1 and evaluating the higher powers of J gives,, 



CO 4 -Jll LQ^ -2&LU Z + JI6U3 4- (oON - 0 . (II-A-2 . 1) 

By requiring that both the real and imaginary parts of Equation II-A- 
2.1 go to zero independently, 

4 2. 

60 - 2660 + 60 N = 0 (II-A-2. 2} 

and 



— 1 1 us* + j 1 6 LO = o 



(II-A-2. 3) 



Dividing Equation II-A-2. 3 by J (JO , gives 

11 UJ 2 - - 16 = 0 



(II-A-2. 4) 



from which the angular frequency of the limit cycle is obtained as 

LL) = 1.206 rad/sec. 

2 . 

Substitution of 60 from Equation II-A-2, 4 into Equation II-A-2. 2 

and solving for N yields 
ave J 

N = 0.595 ♦ 
ave 

To determine the approximate amplitude of the limit cycle, Equations 

II-A-1.11 through II-A-1.15 are again used. Thus, 

£ = 12.402 Volts 

max 

A graphical presentation of the root locus solution is given in 
Figure II-A-2. 1. 

The root locus method predicts an unstable limit cycle by the follow- 
ing reasoning. If the input signal is small, N is small and the loop gain 
is small placing the instantaneous root at ^ ^ in Figure II-A-2. 1. The 
system is stable and the oscillations decrease in magnitude. As the 
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Figure II-A-2.1 

Root Locus Solution for Dead Zone in Error Channel, 

Loop Transfer Function with No Zeros 

oscillations decrease, N also decreases, driving the root location farther 
into the left half plane. Consequently, the system is always stable and 
the oscillations eventually die out. If the input signal is large, N is 
also relatively large making the loop gain high. The instantaneous root 
location is at The system is unstable and the oscillations increase. 

As the oscillations increase, N increases toward N = 1. The root location 
moves farther into the right half plane.. Thus the system is always un- 

i 

stable and the oscillations will increase indefinitely. 
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II-A-3. Describing Function Method. 



The describing function for a dead zone nonlinearity is 



Go — "Sin 'R - R)/ l-R 2 ]] /0_ (II-A-3.1) 

with 

R - (IX-A-3.2) 

where b is the total minus to plus dead zone* 

The linear loop transfer function for the given system in frequency 
response form is 



G (ooj) — 



3.7r 

Jua(JU) + lX0,5jCLH./J(0,/25'JLU+l) 



(II-A-3. 3) 



Table II-A-3.1 gives the magnitude and phase of Equation II-A-3. 3 
for the different frequencies which were used to plot the linear loop 
transfer function on the gain-phase plane of Figure II-A-3.1. These 
values are obtained from a Bode diagram plot of Equation XI-A-3. 3. 



TABLE II-A-3.1 



(JO 


| GOtO) d b 


7C?(J (AJ) cieq 


170 


7.7 


-168 


1.1 


6.2 


-174 


1.2 


4.7 


-179.5 


1.3 


3.5 


-134 


1.4 


2.0 


-189.5 



TABLE II- A- 3. 2 



R 


Go 


*/ 6d 


^ /G 0 db 


0 


1 


1 


0 


0.2 


0.746 


1.34 


2.55 


0.3 


0.624 


1.60 


4.10 


0.32 


0.61 


1.64 


4.29 


0.35 


0.577 


1.73 


4. 73 


0.4 


0.504 


1.98 


5.94 



Table II-A-3.2 gives the values of 1 G/D, in the area of interest. 

d 

These values were obtained from equation II-A-3.1. 
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By interpolation of Figure II-A-3.1, the intersection and thus the 



limit cycle occurs at 



and 



(jJ = 1.21 rad/sec . 



R = 0.34 # 

From Equation II-A-3.2, the amplitude of the limit cycle is then 

£• max ~ 2R 



or 

= 14.71 Volts 

max 

The describing function method predicts an unstable limit cycle by 

the following reasoning. For small input signals, R is relatively large. 

Assuming the -1/G^ curve to be the locus of the system critical pointy then 

the instantaneous critical point is at R in Figure II-A-3.1 and the sys- 

a 

tern is stable due to a positive phase margin. The oscillations will de- 
crease causing R to decrease, thus moving the critical point farther up 
the -180° axis. This motion of the critical point increases the positive 
phase margin. Thus the system is always stable for small signals and all 
oscillations will eventually die out. For large input signals R is small 
placing the instantaneous critical point at R^. The system is unstable due 
to a negative phase margin and the oscillations increase. As the oscilla- 
tions increase R becomes even smaller moving the critical point down the 
-180° axis and increasing the negative phase margin. Consequently, for 
large signals the system is always unstable and the oscillations increase 
indefinitely. 
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II-A-4 Analog Computer Simulation. 

The analog simulation of the given nonlinear system is shown in 
Figure II-A-1. The coefficient pot settings with associated resistances 
and capacitances for real time and real magnitude scaling are given in 
Table II-A-4. 1. 



TABLE II-A-4. 1 



Pot Setting 


Associated Elements 


E. 




a i = Too 


For Input Step of E, Volts 




Rj = 1 Meg, R fl = 1 Meg 


a 2 = 1.0 


R 2 = 1 Meg, R fl = 1 Meg 


= 1.0 


R 3 = 1 Meg, R f2 = 1 Meg 


3 4 = 1-0 


= 1 Meg, = 1 uf 


a,. = 1,0 


Rg. = 1 Meg, C-. = 1 uf 


a^ = 0.6 


R & = 0.1 Meg, C f2 = 1 uf 


a ? = 1.0 


= 0.5 Meg, C ^ 2 - 1 uf 


o 

*— H 
11 

CO 


R g = 0.1 Meg, C = 1 uf 


a 9 = 0.8 


R 9 =0.1 Meg, C f3 = 1 uf 


a io ■ L0 


R 10 = 1 Me§ ’ C f4 = 1 uf 


a ll ■ uo 


R 11 = * R f 3 = * 


a = 0. 5 


For a 5 Volt Dead Zone Each 


Initial Setting 


Side of Zero 



To accurately set the back-bias pots "a n s a sine wave generator was 
connected to the input of the dead zone simulator. An X-Y recorder was 
then calibrated for 5 volts per inch on each axis. The input and output 
of the simulator was then connected to the X-Y recorder and pots were 
adjusted until both break voltages of the static characteristic curve were 
5 volts. The static characteristic curve appears as Figure XX-A-4.2. 
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Ana log Computer Simulation o£ Dead Zone, in the Error Channel s Loop Tcsns£er Function wi 



The system was first excited with a 15 volt step input signal and 
the system was found to be stable as shown in recorder trace II-A-4,1. 

The system was then excited with a 17 volt step input and was found to be 
unstable as shown in recorder trace II-A-4.2. The unstable limit cycle 
was found to occur for an input step of 17o07 and is shown in recorder 
trace II-A-4.3. From recorder trace II-A-4.3 the limit cycle is seen to 
have an angular frequency of 1.207 rad/sec and an amplitude of 15 volts. 
The phase protrait of the system is shown in Figure XI-A-4.3. 



53 




54 





}N CLEVELAND, OHIO PRINTED IN U S A 



Recorder Trace TI-A-4", 1 

Dead Zone in Error Channel, Loop Transfer Function with No Zeros, 15 Volt 

Step Input 

SCALE: 0-1 Volt/Line - 1 Vo It /Line X - 1 Volt/Line 

Tlnie - 1 Second/Divisicn 
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( 





Recorder Trace II-a-4,2 



% 



Dead Zone in Error Channels, Loop Transfer Function with 

No Zeros, 17 Volt Step Input 

SCALE : £-2 Volts/Line £. — 2 Volts/Line * 

X— 2 . Veits /Li ns .Time - 1 Second /Division: 
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mTHiil illliiH 




/ 







Limit Cycle*, Dead Zone in Error Channel „ Loop Transfer Function 
with No Zeros, 16,07 Volt Step Input 

SCALES : 

0-1 Volt /Line 0 - 1 Volt /Line X - 1 Volt /Line 

Time - 1 Second/Division 
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I 



j 

\ 



a' 




H i Ptf |ASc PfiRTRMT - | D?A q -ZONE IN ERROR CHANNEL, L .^_ .60 

LOOP TRANSFER FUNCTION WITH NO ZEROS : ^ S(S+lYS+2)(S+8) 

i l __J ! i j ; 
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II-A-5 Comparison of the four analytic methods in limit cycle prediction 
for the system of Figure II-A # 



Method 


Frequency 

(rad/sec) 


Amplitude 

(volts) 


Mitrovic 


1.206 


12.402 


Root Locus 


1.206 


12.402 


Describing Funct ion 


1.21 


14.71 


Analog Simulation 


1.207 


15.0 
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II-B Loop Transfer Function with One Zero. 

The block diagram of the system chosen to analyze is shown in Fig. 

II-B. 




Figure II-B 

Dead Zone in Error Channel Loop 
Transfer Function with One Zero 
II-B-1. Mitrovic's Method. , 

If N is the instantaneous gain of the nonlinearity the character- 
istic equation for the system is , 

S r tl/.rs\3l.S'$ 3 + : *</s i + (8+5’A/)s + 2 0H=O . (ii-b-1.1) 

Assuming the last two coefficients to be variable, the characteris- 
tic equation becomes, 

S r * II. rs 4 + 3l.sTs 3 + 2^ s 1 + 6 ) S46o-0 (ii-b-i.2) 

where 

i 

8, = 8 + (i’i-b-1. 3) 

and 

6c = 20 N . (II-B-1.4) 
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From Appendix A s Mitrovic c s equations for and 

6 0 = -[2 £ / Wft a ^,^)43 ISm0 a (!S) +i \Su n %08) 4 oo/fadf)] (ii-b-i.5) 

and 

6,= 2.9u,^(»)+3/,r Wr ,^3) + l|.i'W B 3 ^)+U)„VfW3 « (II-B-1.6) 

Substitution of the values of the 0 functions from Appendix A for 

%-€> gives, 

6 0 = 29 W„ 2 - I hS-UJ„ 4 (11-1-1.7) 

and 

6, = 3J.rw, l -iO, 4 . (II-E-1.3) 

The sets of parametric equations IX-B-1.3, IX-B-1.4, XX-B-1.7 9 ^nd 
II-B-1.8 describe the M-point locus and the stability curve on the 
versus B. plane. Direct solution of these equations for the intersection 
becomes fairly complicated; consequent ly 9 a straight graphical solution is 

performed. 

Table IX-B-1.1 gives the values of B q and for values of (j 
plotted in Figure XX-B-1.1. Elimination of N from Equations XX-B-1.3 and 
IX-B-1.4 yields 

8i = 8 + -§£- (n-E-1.9) 

4 

which is a straight line with a B.. intercept of 4 8 and a slope of -f 4. 

The intersection of the two curves is seen to be approximately at £) = 0.56 

rad/sec . 

To accurately determine the intersection more values of stability curve 
are calculated and the value of B, thus determined is substituted in Equa- 
tion XI-B-1.9. 
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TABLE II-B-1.1 



U) a 


B 

o 


B 1 


0.5 


6.53 


7.81 


0.6 


8.95 


11.21 


0.8 


13.85 


19.75 


1.0 


17.5 


30.50 


1.2 


17.9 


43.23 


1.4 


12.6 


57.86 


1.588 


0 


73.0 



The intersection occurs when the two values of B q are equal. These addi- 
tional calculations appear in Table II-B~1,2 



TABLE II-B-1.2 





B 1 


B 

o 


B 1 

o 


0.55 


9.43724 


8.72017 


5.74897 


0.567 


10.02355 


8.13459 


8.09419 


0.568 


10.05857 


8.15910 


8.23428 


0.57 


10.12679 


8.20816 


8.50716 



By interpolation the intersection occurs at 

(A)*) = 0,5673 rad/sec 

and 

B =8.140 
o 

From Equation II-B-1.4, the average nonlinear gain over a cycle is 

N = 0.407 
ave 
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Since the nonlinearity is the same as the one of Section A-l of 
this chapter, Equations II-A-1.11 through II-A-1.15 are used to obtain 
the approximate amplitude of the limit cycle. Substituting into 

Equation II-A-1.14 yields 

A = Tan' 1 0.407 = 22.15° . 

Substituting A into Equation II-A-1.15 yields 

B *= 22.85° , 

Substituting A and B in Equation II -A-l. 13 gives 



r sm 22.15 , or 

a = 5 — y- - — = 4.85 , 

sin 22,85 

Substituting "a" into Equation II-A-1.12 gives 

c = 0,707 (4.85) - 3.42 # 
Finally Equation II-A-1.11 gives 



£ = 5 + 3.42 = 8.42 volts # 

max 

Applying the same reasoning as given in Section II-A-i to Figure 
II-B-1.2, Mitrovic's method predicts an unstable limit cycle. 



63 



II-B-2. Root Locus Method. 



Substitution of into Equation II-B-1.1 and simplifying the 

higher order J terms gives 

JCO Jr f/;.5'6o' , -J3l.ra3 3 -2 £ ?^ 2 +J^ + 5'N)a)4iON-C , (II-B-2.1; 

By requiring that both the real and imaginary parts of Equation 

11-6-2.1 go to zero independently gives 

oc/-j3i,^o*V,)(&+,5>i)u) - o Co-2-2 .2) 

and 

\1>S UJ“ 2.9 LQ+2CW - 0 ♦ (II £-2.3) 



Dividing Equation II-B-2 C 2 by Ju) and solving for 5h gives 

5 N = - + 31 . 5 tO -8 » {II-B-2.4) 

Substitution of Equation XI-B-2.4 into Equation H-E-2.3 and simplify- 
ing gives 

7.5 (JO 4 + 97 UJ 2 - 32 = 0 , 

2 

Solving this equation for UJ gives 



Uj = 0.3219 



from which 



UJ = 0.567 rad/ sec . 

2> 

Substituting the value of MJ ^ back into Equation II-E-2.4 gives the 

’’average* 1 nonlinear gain over a cycle as v 

N - 0.4072 . 

ave 

To find the approximate amplitude of the limit cycle ^ Equations 

II-A-1.11 through IX-A-1.15 are again used. These equations yield 



£ - 8.41 volts 

max 
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2 . 1 . 



A graphical presentation of the solution is sketched in Figure II-B- 



Applying the same reasoning as given in Section II-A-2 to Figure II- 



B-2.1, the root locus method predicts an unstable limit cycle. 




T T T, _ 0 T* 

r i i - c - 



Root Locus Solution for Lead Zona In t n z- 



»- /'■V . 



Loop Transfer Function with otd Z- o 
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II-B-3. Describing Function Method, 

The describing function for dead zone is given in Equation XI-A-3*!. 
The linear loop transfer function for the given system in frequency 

response form is 



/- / . „ 2.r(o.2rotjj -h) . 



(II-B-3. 1) 



Table II-B-3. 1 gives the magnitude and phase of Equation II-B-3. 1 
for the different frequencies which were used to plot the linear loop 
transfer function on the gain-phase plane of Figure II-B-3. 1. These 
values are obtained from the Bode diagram plot of Equation IT-E-3.1. 



TABLE II-B-3. 1 



CO 


) Goto) | db 


/GPfrOUeq 


0.495 


+ 10 


-172 


0.53 


+ 9 


-176 


0.56 


+ 8 


-130 


0.60 


+ 7 


-134 


0.638 


+ 6 


- 138 



Table II-B-3. 2 gives the values of i/G plotted in Figure II-B-3. 1. 

u 

These values are obtained from Equation II-A-3.1 for various values of R„ 
, TABLE TX-B- 3. 2 



R 


Go 


1/Gd 


J/Go ctb 


0.4 


0.5044 


1.983 


5.94 


0.45 


0.447 


2.237 


6.99 


0.48 


0.4132 


2.420 


7.48 


0.49 


0.4021 


2.4817 


7.91 


0.5 


0.3910 


2.558 


3.16 



By interpolation of Figure II-B-3. 1 the intersection and thus the 
limit cycle occurs at 
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CO -0.56 rad/ sec 



and 

R - 0.495 ♦ 

From Equation IX-A-3.2, the amplitude of the limit cycle is then 



£ - 

max 


2(0.495) * “- 1 volts ' 



Applying the same reasoning given in Section II-A-3 to Figure II-B 
3.1, the describing function predicts an unstable limit cycle. 
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XX-B-4 Analog Computer Simulation. 



The analog simulation of the given nonlinear system is shown in 
Figure XX-A-4.1. The coefficient pot settings with associated resist 
ances and capacitances for real time and real magnitude scaling are g 

in Table IX-A-4.1. 



TABLE XX- A- 4.1 



Pot Setting 


Associated Elements 


E „ 

1 

*1 “ 100 


For Input Step of volts 

R„ = 1 Meg lS , R fl = I Meg 


a 2 = 1.0 


1>2 8 1 Meg, R £| S 1 Mfcg 


a 3 = 1.0 


R 3 = 1 Meg s R f2 = 1 Meg 


a. = 0.4 
4 


R, - 0.4 Meg, €*,„ - 1 of 


= 0.8 


= 0.1 Meg, “ 1 uf 


a = 1.0 
6 


Rg “ 1 Meg, R f3 - 1 Meg 


a = 1.0 
7 


R - 1 Meg, ~ 1 Meg 


a 8 = 0.5 


Rg = 0.1 Meg, C f2 » 1 uf 


a 9 = 1 *° 


R 9 *= 0.5 Megs, C f2 « 1 uf 


a iA “ 1>0 

10 


R 10 = 1 M ' ga C f3 = 1 uf 


a s 1.0 

11 


R n = 1 Meg, C f3 - 1 uf 


a^2 ~ 1.0 


R ^2 “ 1 Megj, = 1 af 


*13 " °' 5 


Rjg = 1 Meg, C f4 s 1 uf 


a X4 “ ^ 0 ® 


R u = 1 Meg, C f5 = 1 u£ 


a = 0.05 
Initial setting 


For a 5 volt Dead Zone 
Each Side of Zero 


The back-bias pots ! 


! a ?} were accurately set " in the same manner 



described in Section IX-A-4. The static characteristic of the non- 
linearity is shown in Figure XX-A-4.I. 



.van 
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+100 v» 
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The system was first excited with a 10.2 volt step input and was 
found to be stable as shown in recorder traace II*B-4.1. The system 
was then excited with a 10.6 volt step input and found to be unstable* 
as shown in recorder trace II-B-4.2. The unstable limit cycle was 
difficult to obtain but was found to lie between step inputs of 10.3 volts 
and 10.4 volts. These two responses are shown, in recorder traces TI-£-4.3 
and XX-R-4.4. The average of the two traces show the angular frequency to 
be 0.552 rad/sec and the amplitude*, using the first two peaks of each 
trace* to be 9.9 volts. 

The phase portrait of the system is shown in Figure TI-B~4«2. 
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Recorder Trace II-B-4.1 

Dead Zone in Error Channel, 
Loop Transfer Function with 
One Zero, 10.2 Volt .Step 
Input. i 

SCALES: 

- 1 Volt/Line 
C. ** 0.5 Volts/Line 
Time - 5 Seconds/Division 




ISTRUMENTS division c 




./ 



a 



Recorder Trace II-B-4.2 

Dead Zone in Error Channel, 
Loop Transfer Function with 
One 7ero, 10.6 Volt Step Input 

SCALES : 

£ - 1 Volt/Line 
&. - 0.5 Volts/Line 
Time - 5 Seconds/Division 
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Recorder Trace XX-E-4.3 

LIMIT CYCLE, Dead Zone in Error Channel, Loop Transfer Function with One Zero 
SCALES: £ - 1 volt/line X - 0*5 Volt/lime Time - 1 Second/division 
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II-B-5 Comparison of the Four Analytic Methods in Predicting Limit 
Cycles for the System Shown in Figure II-B. 



Method 


Angular Frequency 
rad/ sec 


Amplitude 1 

Volts 1 


Mitrovtc 1 s 


0.5673 


8.42 1 


Root Locus 


0.567 


8.41 1 


Describing Function 


0.560 


10.1 1 


Analog Simulation 


( 0.552 


9.9 1 
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CHAPTER III 

IDEAL RELAY IN ERROR CHANNEL 



I 



III-A Loop Transfer Function With No Zeros 

The system arbitrarily chosen to analyze is show in the block 
% v 

diagram of Figure III L A. 



/ i 



i 




Figure III-A 

Ideal Relay in Error Channel 
Loop Transfer Function With No Zeros 



III-A- 1 Mitrovic's Method 

Assigning N as the instantaneous gain of the relay, the character- 
istic equation can be written as 



S^tV.Ss 4 +?>0.£s 3 t-375 2 + I2.S -H2N =0 . 



(III-A-1.1) 



Assuming the last two coefficients variable, 



the 



characteristic • 



equation becomes, 

30.5'£ 3 + 37s 2 ' 4 8 , 5 + Bo - (III-A-1.2) 

where 



B = 12 N 

o 



(III-A-1. 3) 



I 
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and 

B. = 12 0 (m-A-1.4) 

X 

From Appendix A s Mxtrovic rj s equation for B and a ra $ 

O A. 

Bo - -[37(o n a (^£<j) +3 i LdffaCi)] (iii-a-i.3> 

and 

6, = 3 7 W „ ^ l^ + SOSto^ «0 +9, . (II1“A-1.6> 

Substitution of the values of the 0 functions from Appendix A for 
^ s o gives the parametric equations of the stability curve as s 

B o = 37 &V, 2 - 9.5 UJ n 4 (IlI“A“1.7ii 

and 

B x = 30.5 u) n 2 - 60^ 4 t (III-A-1.8) 

To find the intersection between the M- point locus (Equations XXX-A- 
1.3 and XIX-A-1.4) and the stability curve (Equations XXI-A-l*? and XIX- 
A-1.8) s the graphical approach is chosen. Table ITI-A-l.l gives the 
values of and for various U0 n that are plotted in Figure TXI-A-l o l 0 
The M-point locus is seen to be a straight line parallel to the axis 
at - 12. From Figure XXX-A-i.i y the intersection ’is seen to occur at 
about 60^- 0.62. 
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TABLE III-A-1.1 



60 k) 


B 

o 


B 1 


0.5 


8.65 


7,56 


0.6 


12.08 


10.85 


0.7 


15.85 


15.70 


0.8 


19.8 


19.2 


1.0 


27.5 


29.5 


1.2 


33.6 


41.8 


1.4 


36.0 


56.0 


1.6 


32.5 


71.4 


1.8 


20,3 


88.2 


1.974 


0 


103.6 



To accurately determine the intersection more values of B and 

o 

B 1 are calculated for U>t\ close to 0.62. These values appear in Table 
III-A-1.2. By interpolation from Table III-A-1.2 the intersection and 
thus the limit cycle occurs at 

^■0^ = 0.63X5 rad/sec 

and 



B = 13.24 
o 



TABLE III-A-1.2 





B 

o 


B 1 


0.62 


12.8191 


11.5764 


0.63 


13.1891 


11.9480 


0.631 


13.2259 


11,9854 


0.632 


13.2631 


12.0229 



The "average" nonlinear gain over a cycle is determined from 
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Equation III-A-1.3 as > 



N - 1.103 



ave 

From this the approximate amplitude of the limit cycle is 

= 4.53 volts 

max 

Note that the M-point of the linear system in Figure ITX-A-l.l shows 
that the linear system is stable. It is well known that am ideal relay 
gives a limit cycle for a linear system of order three or above and the 
occurance of a limit cycle is independent of linear system stability. 

This is* easily explainable on the versus plane. The ideal relay 
has an infinite maximum instantaneous gain whereas the nonlinearity used 
in Chapters I and II have a maximum instantaneous gain of one. Thus on 
the versus B^ plane the M-point locus is unlimited for the ideal re- 
lay but does have a maximum value for the other nonlinearities. Conse- 
quently, an infinite straight line segment must necessarily intersect the 
stability curve at some point producing a limit cycle. 

Mitrovic's method predicts a stable limit by the same reasoning used 
in Chapter I, Section I-A-l, This same argument can be applied to the 
ideal relay since this nonlinearity can be envisioned as saturation with 
an infinite gain in the linear region. 
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XII-A-2 Root Locus Methods 

Substituting S =Ju) into the characteristic Equation ‘fir -A-?. I a:.i 
evaluating the higher older powers of J yields,, 



^+^60 4 -J36.5"60 3 -37^-f j iaoJ + /aM = C, 



Requiring that the real and imaginary parts of Equal ion IIX-A-l*! 
go to aero independent ly gives, 

Jui' s '-j30.vT6 jo 3 -tW-u. j-0 t/rx-A-2 

and 

9.v5’u) 4 -37o> x -f i C N > {nr-A-2 

Dividing Equation XII-A-2.2 by and .solving for U) gives 



from which 



do = 0,3986, 



d>0 - 0,6314 rad/sec ^ 



Substituting 00 into Equation 11I-A-2 0 3 and solving for h gives, 

N - 1 c 1032 
ave 

From this "average** nonlinearity gain over a cycle, the approximate 

amplitude is determined as, 

£ = 4. 533 voles 

max 

The graphical presentation of the solution by root 1 ecus .methods 

is shown in Figure XIX-A-2,1, 

The root locus method predicts a stable Unit cycle fey appl/i g -he 
same arguments to Figure III-A-2,1 that are given in Chapter lh, Section 

I-A-2. 
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/ 




Fig* 

Root Locus Solution for Ideal j&ei&y in L-ror ufocnnei 

i 

Loop Transfer Function with No Zeros 

t 



8 d 



IIT-A-3 Describing Function Method. 



The describing function for an ideal relay is 

K o 



g„ - 

D ITS 



LSL 



(llI-A-3.1; 






The linear loop transfer function for the given system in frequency 
response form is, 



0W(2.J(^+|)(O,£jUJi-l)(O.333^UJ+/XO'^^+0 » (IXT-A-3. 2) 

Table IXI-A-3,1 gives the magnitude and phase of Equation 11I-A-3.2 
for the different frequencies which were used to plot the linear loop 
transfer function on the gain-phase plane of Figure IXI-A-3.1* These 
values are obtained from a Bode diagram plot of Equation IXT-A-3.2. 



TABLE IXI-A-3.1 



(JO 


| G0w)| ^ 




0.527 


2.0 


- 167.5 


0.565 


1.0 


-172 


0.605 


0 


•176.5 


0.625 


-0.5 


-178 . 5 


0.648 


-1.0 


“181 


0.69 


-2.0 


•185.5 


0.735 


-3.0 


-190 



Table IXX-A-3.2 gives the values of 1/G in the area of interest as 
obtained from Equation III-A-3.1 for various values of error. 

TABLE XXI-A- 3. 2 



£ rAifl y 


&n 


l/Gn 


j/ G jo db , 




1.273 


0.7854 


- 2.10 


6 


1.061 


0.942 


-0.52 


5.8 


1.098 


0.911 


-0.80 



34 





By interpolation of Figure IXX-A-3.1 S the intersection and thus 



the limit cycle occurs at 

00 - 0.635 rad/sec 

and 

C. = 5«b volts o 
max 

The describing function method predicts a stable limit cycle by 
applying the same arguments in Section X-A-3 of Chapter I to Figure XIX- 
A-3.1 
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III-A-4 Analog Computer Simulation. 

The analog simulation of the given nonlinear sy&te.r is shews in 
Figure III-A-4. 1. The coefficient pot settings with associated resist- 
ances and capacitances for real time and real magnitude scaling are given 
in Table III-A-4. 1. 



TABLE III-A-4. 1 



Pot Setting 


Associated Elements 


E , 




a = — — 


For Input Step of S. Voir. 


1 100 


Rj - 1 Meg, R fl = 1 1 Meg 


a ? - 1.0 


R 2 = 1 Meg, R fl = i Meg 


a^ ~ (See discussion) 


R.^ = 0.1 Meg, R f2 = 1 Meg 


a 4 = 0.4 


R 4 = 0.1 Meg, C f1 » 1 u£ 


a„ - 0.5 
0 


R s » 1 Meg, C fl - 1 uf 


a 6 = °” 3 


R, "0.1 Meg, C = 1 uf 

o 


a_ = 1.0 
7 


= 0.5 Meg, ~ 1 uf 


a 8 = 1,0 


R g = 1 Keg, C f3 « 1 uf 


A g ~ 0.3 


R g , = 0.1 Meg, C f3 = 1 uf 


&> 

o 

tt 

►-* 

o 


^10 ~ 3 hieg, C^ 4 - 1 uf 


^11 ~ ® ^ 


R n - 0.1 Meg, C f4 = 1 uf 


a l2 s 1 


Rj 2 = 1 Mag, C fe . ~ 1 uf 


a = 0.05 (Initial 


For 5 Volt Voltage 


sett ing) 
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Figure. III-A-4.1 

Analog Computer Simulation of Ideal Relay in Error Channel. Loop Transfer Function with No Eeros. 



The simulation of the ideal relay is seen to be a simple saturation 



simulation. However* the gain of the linear portion is made very high 
by the insertion of coefficient pot a^ 0 To set the back-bias pots 
and pot & sine wave generator was connected to the input of the re- 
lay simulator. The output was monitored and the back-bias pots were 
then adjusted until the output was limited at ± 5 volts. The output and 
input was then connected to an X - Y recorder to obtain a static character 
istic curve. Pot was then adjusted until the slope of the linear por- 
tion was as nearly vertical as possible without the amplifier overloading. 
The static characteristic curve appears as Figure IXX-A-4.2. 



The system was first excited with a 1 volt input signal and was found 
to build up to the limit cycle as shown in recorder trace XXX-A-4.1. The 
system was then excited with a 10 volt step input and was found to go into 
the same limit cycle as shown in recorder trace XXX-A-4.2. From these two 
recorder traces the limit cycle is seen to be stable with an angular fre- 
quency of 0.628 rad/ sec and an amplitude of 5,7 volts. 

The phase portrait of the system is shown in Figure IXI-A-4o3. 
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Recorder Trace IXX-A-4.1 



Ideal Relay in Error Channel , Loop Transfer Function with No Zeros 
1 Volt Step Input 

SCALES ; 

- 0o5 Volts/ Line 



8- - G. 5 Volts/Line X - 0.5 Volts/lint 

Time * 5 Seconds/Division 
1 Second/Divi sicn 
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\ 






f : 




SCALES: 

8 - 0,5 Vo Its /Line 



# 

8 - 0,5 Volts/Line x 

Tinie - 5 Seconds /Division 
1 Secoind/Divisicn 



,5 Vol' 
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IXX-A-5 Comparison of the Four Methods in Predicting Limit Cycles for 
the System of Figure IIX-A. 



Method 


Angular Frequency 
rad/ sec 


Amplitude 

Volts 


Mitrovic * s 


0.6315 


4.53 


Root Locus 


0.6314 


4.533 


Describing Functions 


0.635 


Kf Q 
. J * c 


Analog Simulation 


0.628 


5 - 7 
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III-B Loop Transfer Function with One Zero. 

The system to analyze is shown in the block diagram of Figure 
III-B. 

» 

•< T 




Figure III-B 

Ideal Relay in Error Channel Loop 
Transfer Function with One Zero 
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III-B-1. Mitrovic's Method. 



Assigning N as the instantaneous gain of the relay, the character- 
istic Equation can be written as 

5 6 -f -31 s 5 1 3IOs 4 + 1^40 S 3 f m<\ sN (1024 +$ooon)s i4000N=O. 

(IXI-B-1.1) 

Assuming the last two coefficients are variable* the characteristic 
equation becomes* 

S^+ 3/S^+3l0s 4 +i2405 3 -f I984s a -+3| 5 4 6o = 0 



where 



and 



B --4000 N 
o 






Bj = 1024 4- 8000 N (XII-B-1,4) 

From Appendix A* Mitrovic°s equations for and B 1( are., 

8 0 = -[mm uv\0,(g)+ i2Ww,^j(if)-i 3l6io„ 4 </ 3 (i) 

+ 3liO^(H)4 U0„ k </f(y)J nn-B-1.5) 

and 

8, = I9?4 0j„4i(/-i) + t 3«to,M W) 

+ 3lu)„%-(«)+W«' r <&(S0 - 

Substituting the values of the functions from Appendix. A for ^ - 0 



gives the parametric equations of the stability curve as* 



8 0 = /984 



(111-E-l . ?) 



and 

8, - 13^0 6o n x - 3/ o)ki 4 



(III-B-1.3) 



To find the intersection between the M- point laiue (Equations IYJ-B-I.3 
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and II1-B-1.4) and the stability curve,, the graphical approach is chosen. 

Table III-B-1,1 gives the values of B and B, for the stability curve as 

o 1 

plotted in Figure III-B-1. 1. 



TABLE III- B- 1.1 



Uj ^ 


B 

0 


B i 


0.5 


227 


308 


0.8 


1143 


7&0 


1.0 


1675 


1209 


1.2 


2219 


1722 


1.4 


2708 


2311 


1.6 


3067 


2962 


1.8 


3200 


3695 


2.0 


3040 


4464 


2.2 


2454 


5275 


2.3 


1968 


5692 


2.5 


544 


6540 



To obtain the M- point locus solve Equation ITT-B-1.4 for N and 
substitute into Equation 1II-B-1.3. 

B o = °* 5 B x “ 512 (III-B-1 . 9) 

Equation III-B-1. 9 is seen to be a straight line with a slope of 0.5 and 
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a B, intercept of 1024. 

It Is seen from Figure III-B-1.1 that the Intersection occurs at 
roughly = 2.24. To more accurately determine the intersection, 

values of B and B, are calculated for LG i a close to 2.24. The value 

o 1 n 

of B. is then substituted into Equation III-B-1.9 and the intersection 
1 

occurs when the two B's are equal. These additional values are given 
in Table III-B-1.2. 



TABLE III-B-1.2 



i On 


B (Stability 
curve) 


B 1 


B (M- point 
Locus) 


2,24 


2266 


5439 


2208 


2.249 


2233.5 


5478.8 


2227.4 


2.25 


2228.7 


5483.0 


2229.5 



Thus from Table III-B-1.2 # the frequency of the limit cycle is 

£d n = 2,25 rad/sec ? 

Substituting B q at = 2.25 into Equation III-B-1.2;, the "average" 

gain over a cycle is determined as 



N 



ave 



2229 

4000 



from which the approximate amplitude is 



= 0.5573 , 
determined as 



£ 



max 



5 

0.5573 



8.973 volts. 
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Mltrovic's method predicts a stable limit cycle by the same argu- 
ments applied to Figure III-B-1.1 as were given in Chapter I„ Section 
I-A-l. 
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_L 
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III-B-2 Root Locus Method. 



Substituting s = j(jJ into the charaecerist ic Equation III-B-l.J 
and evaluating the higher powers of J yields* 

31 to 5 " + 31 0 U4 4 -J alOud 3 -i9?4w ^ 

+Xl&M + g060M)lO+4000KJ — O. 'm-S-2.1) 

Requiring that the real and imaginary parts of Equation 111*2.-2.1 

go to zero independently* gives 

-do L + 3 1 0 to 4 - I « £4 £o 2 + 4060 N = 0 O it. - b- 4 * , 

and 

J3I oc^-U 1*464 3 + ^0^04 + 8^00V)i0-O * (in-a-i.r. 

Dividing Equation III-B-2 ,3 by and solving this for LQOQ K 

gives* 

4A30N - 620 w 2 ->T/a- if.rto 4 - (tii-b-:.*.,. 

Combining Equations III-B-2. 3 and III^B-2.4 gixes 

<0^-2Hr’6o 4 + /3&to% r/1 -o. (ITI-5-2. 5) 

2 

To solve for 6-0 from Equation III L B-2.5;, a digital computer pro- 
gram was used and gave the roots> 

Ui 2 = 5.062 

(O 2 = - 0.55 

<*t) 

(JO^ -299.55 

Only the first value is of use and gives 

(jj “ 2.50 rad/ sec. 

2 

Substituting (jO into equation III-B-2. 4 gives the 'average” gain 

of the relay over a cycle as 

N *= 0.5574 
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from which the approximate amplitude of the limit cycle is determined to 



£ = 8.972 volts 

max 

A graphical presentation of the solution is sketched in Figure 

% 

rtl-B-2.1. The root' locus method predicts a stable ( limit cycle by 
applying to Figure III-B-2.1 the same arguments given in Chapter I, 
Section II-A-2. 




Figure III-B-2.1 

Root Locus Solution for Ideal Relay in Error Channel, 
Loop Transfer Function with One Zero. 
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IXI-B-3 Describing Function Method, 

The describing function for the ideal relay is given in Equation 
TXI-A-3,1, The linear loop transfer function for the given system in 
frequency response form is 

r/)M)- , ,, 3.9/(2jtUfQ 

Table IXI-B-3.1 gives the magnitude and phase of Equation III-B~3.I 
for the different frequencies which were used to plot the linear loop 
transfer function on the gain-phase curve of Figure III-B-3.1. These 
values are. obtained from a Bode diagram plot of Equation TII-E-3.1. 



TABLE III-B-3.1 



to 


|GQto)| <±k 


/ ciee 


1.97 


7 


-165 


2.1 


6 


-174.5 


2.24 


3 


-180 


2.38 


4 


-185 


2.52 


3 


-190 



Table III-B-3.2 gives the values of 1/Gj, in the area of interest as 
obtained from Equation XXI-A-3.1 for various values of error. 

TABLE XIX-B-3.2 



E 

max 


S 


I/e D 


i/Mb 

L> 


8 


0.796 


1.257 


1 . 99 




10 


0.634 


1.571 


3.93 




11 


0. 579 


1.725 


4.75 




11.3 


0.563 


1.775 


4.98 




12 


0.531 


1.885 


5.51 
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By interpolation of Figure III-B-O.l*, the intersection and thus 



the limit cycle occurs at 

6 0 - 2,24 rad/sec, 

and 



E a 11,31 volts # 
max 

The describing function method predicts a stable limit cycle by 
applying to Figure IIX-B-3.1 the same arguments given in Chapter I* 
Section I-A-3. 
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XII-B-4 Analog Computer Simulation* 

' The analog simulation of the given nonlinear system is shown in 
Figure IXX-B-4.1. The coefficient pot settings with associated re- 
sistances and capacitances for real time and real magnitude scaling are 
given in Table IXI-B-4,1. 



TABLE XXX-B-4* 1 * 



Pot Setting 


Associated Elements 


E. 


For a Step Input of E, Volts 


a i ~ 100 


Rj = 1 Meg j, R fl = 1 Meg 


~ 1*0 


R 2 = 1 Meg s R^,, = 1 Meg 


= (See discussion) 


R 3 = 0.1 Meg, R = 1 Meg 


a =0.5 
4 


R 4 «* 1 Meg s G fl = 1 uf 


a 5 = 1.0 


R 5 = 1 Meg s C fl = 1 uf 


a, = 1.0 
6 


Rg = 1 Megs R f3 = 1 Meg 


a^ = 1.0 


R ? = 1 Megs R f3 - 1 Meg 


a g = 1.0 


R g = 0.1 Megs " 1 uf 


a 9 = 1.0 


R 9 = 0.5 Megs C f2 = 1 uf 


a 10 = 1 -° 


R 10 = 0.1 Meg s C f3 = 1 uf 


a n = 0.4 


R u = 0.1 Meg, C f3 - 1 uf 


= 1.0 


R i2 = 0.1 Meg, = 1 uf 


*13 = 0.8 


R 13 = 0.1 Megs C f4 = 1 uf 


*14 = °' 8 


R i4 » 1 Megs C f5 = 0.1 uf 


*15 * °' 8 


R T j, - 0.5 Megs - 0.1 uf 


*16 “ K0 


R„, - 1 Megs C = 1 uf 

16 to 


a 17 = 1.0 


R 17 » 1 Megs R f4 * 1 Meg 


a = 0.05 (Initial 


For 5 Volt relay voltage 


Setting) 
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^ s+lk S+g- ' • S-M 

Figure IXI-B-4,1 

Analog Computer Simulation of Ideal Relay in Error Channel, Loop Transfer Function with One Zero. 



The back-biasing pots "a" and pot were set in the same way 
described in Section XXI-A-4. The static characteristic curve is also 
the same as shown in Figure XXX-A-4.1. 

The system was first excited with a 5 volt step input signal and 
was found to build up to the limit cycle shown in recorder trace IXX- 
B-4.1. The system was then excited with a 20 volt step signal and 
found to go into the same limit cycle as shown in recorder trace IIX-B- 
4,2. From these two recorder traces the limit cycle is seen to be stable 
with angular frequency of 2.22 rad/sec and an amplitude of 11.7 Volts. 

The phase portrait of the system is shown in Figure XXI-B-4.2. 
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Recorder Trace ITI-B-4.2 

Ideal Relay in Error Channel, Loop Transfer Function with One Zero, 20 Volt Step Input. 

SCALE o . £. 2 Volts/Line 1 Volt/Line X ~ 0,5 Volts/Lme Tifue - 1 Second/Division 

' 0,2 Second/Division 
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III-B-5. Comparison of the Four Methods in Predicting Limit Cycles for 
the System Shown in Figure III-B. 



r 

Method 


Angular Frequency 
rad/ sec 


Amplitude 

Volts 


Mitrovic 0 s 


2.25 


3.973 


Root Locus 


2.250 


8.922 


Describing Function 


2.24 


11.31 


Analog Simulation 


2,22 


11.7 
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CHAPTER IV 



RELAY WITH DEAD ZONE IN ERROR CHANNEL 
IV-A. Loop Transfer Function with No Zeros. 

% The block diagram of system chosen ife shown in 'Figure IV-A. 




Figure IV-A t 

Relay with Dead Zone in Error Channel Loop 
Transfer Function with No Zeros 

IV-A- 1 Mitrovic's Method 

The characteristic equation of the given system is, 

+ IS’.5'5^+77. < 5'S" 1 + ISS S 2 + 124 5 Z + ?>Z S + SO N - 0. 

(IV-Arl.l) 



Assuming the last two coefficients variable, the characteristic 
equation becomes 

5 S (iTrs 7 7. 5-s*+/5\rs 3 -M2 4 s 2 - + 8, S + 6 0 =0 (IV-A-1.2) 
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where 



B = 50 N 
o 



(XV-A-1.3) 



and 




From Appendix A, Mitrovic's equations for and are, 

+ iSSco^fac# ) + lo^0 s (<4) J 



(IV- A- 1.4) 



(XV-A-1.5) 



and 



B, - HUtoMty + l&S’Uitf'faUf) +7 ~>.&to n 3 <f><iC4) 

4 ir.5fc> n 4 ^D+ u> n s ^us) . 



(IV-A-1.6) 



Substitution of the values of the 0 functions from Appendix A for 
^ = 0 gives the parametric equations of the stability curve as*, 

B q = 124 u; 0 2 - 77.5 U)„ 4 + W n 6 (IV-A-1.7) 

and 

Bj = 155 (Dr, 2 - 15.5 (JO n 4 „ (IV- A- 1.8) 

Table IV-A-1.1 gives the points of stability curve plotted in Figure 
IV-A-1.1. These points are calculated from Equation XV-A-1.7 and IV-A-1,,8 
for various 60^ . 

From Equations IV-A-1.3 and IV-A-1.4 S the M- point locus is seen to 

be a straight line parallel to the B axis with a maximum of B - 50* 

o o 

From Figure IV-A-1.1, the intersection of the M-point locus with the 
stability curve is seen to occur at roughly 60^= 0.46 rad/sec. To more 
accurately determine the intersection more points in the area of the 
intersection are calculated in Table IV-A-1.2. 
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TABLE IV- A- 1.1 



LO n 


B 1 


B 

o 


0.2 


6.18 


4.76 


0.4 


24.0 


17.9 


0.5 


37.8 


26.3 


0.6 


53.8 


34.6 


0.7 


72.2 


42.4 


0.8 


92.9 


48.1 


0.9 


115.4 


50.1 


1.0 


139.5 


47.5 


1.1 


164.8 


38.3 


1.2 


190.9 


25.5 



TABLE IV-A-1.2 



COn 


B 

o 


B 1 


0.45 


21.94 


30.75 


0.46 


22.78 


31.59 


0.461 


22.86 


32.24 


0.47 


23.62 


33.48 



By Interpolation of Table IV-A-1,2 using Equation IV-A-1.4, the 
intersection and thus the limit cycle occurs at 



tOr\ ~ 0.4607 rad/sec 

and 

B - 22.84 # 

o w 
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From Equation IV-A-1.3, the ’’average'* nonlinear gain over a cycle 
is determined as, 

N = 0.457 

from which the approximate amplitude of the limit cycle is determined to 

be £ 10.95 Volts , 

max 

Mitrovic's method predicts a stable limit cycle by the following 
reasoning. The input signal must be larger than the dead zone for the 
system to respond at all. Once the input signal becomes larger than the 
dead zone, the characteristic of the relay is for practical purposes 
similar to saturation. Thus the same arguments given in Chapter I, 
Section I-A-l can be applied to Figure IV-A-1.1. 
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IV-A-1.2 Root Locus Method. 



Substituting 5 = Jlx) in the characteristic Equation TV-A-1.1 and 
evaluating the higher powers of J gives, 



HSui I SSuS 1 - I24L0 1 + )32.C0 + 6~ON - 0 . 

(IV-A-2.1) 



Requiring that the real and imaginary parts of Equation IV-A-2.1 go 
to zero independently, yields the following two equations. 



tO « 0. 210899 

as the only meaningful root. Thus the frequency of the limit cycle is 

UJ = 0 .4592 rad/sec. 

2 

Substituting CO back into Equation IV- A- 2. 2 and solving for N 

as the ’’average" gain of the nonlinearity, gives 

N = 0.4543 
ave 

from which the approximate amplitude is found to be 

£ = 11.01 volts # 

m 

A graphical presentation of the root locus solution is shown in Figure 

IV-A-2.1. 

The root locus method predicts a stable limit cycle by the following 




(IV-A-2.2) 



and 



^ \f,S~ to* - $ ISSco 3 + J3 ILO-O, 



(IV-A-2.3) 
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arguments. For an input signal larger than the dead zone, the relay 
characteristic is similar to saturation and the arguments given in 
Chapter I, Section I-A-2 can be applied to Figure IV-A-2.1. 



\ 






? 



i 



uneaR SYS7FM closed loop root 



LIMIT CYCLE 
Kie-OASH 




Figure IV-A-2.1. 

Root Locus Solution for Relay with Dead Zone 
in Error Channel, Loop Transfer Function with 
No Zeros 
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IV-A-3. Describing Function Method. 



The describing function for a relay with dead zone is 



G _ 4V 



I - (A~) ICL 



where 



(IV-A-3. 1) 

V is the signal controlled by the relay 
^^is the relay input signal 

is the total dead zone minus to plus. 

The linear loop transfer function for the given system in frequency 
response is 

/, 5 % 2 . 

' ' (IV-A-3.2) 



(2 lj_5^2- 

K J Jw^2Ju)+/jptt)+iXasjwtiXo.zrjto+i)(o,i2rjoi+i) 



Table IV-A-3, 1 gives the magnitude and phase of Equation TV-a- 3.2 for 
the various frequencies that are used to plot the linear loop transfer 
function on the gain-phase plane of Figure XV-A-3.1. These values are 
obtained from a Bode diagram plot of Equation IV-A-3, 2, 



TABLE IV-A-3, 1 



60 


|&0o)| Ub 


/&{M) de<) 


0.4 


9 


-172 


0.43 


8 


-177.5 


0.455 


7.5 


-180 


0.463 


“9 

/ 


-182 


0.49 


6 


-186.5 



Table IV-A-3.2 gives the values of 1/G^ which are plotted in 
Figure IV-A-3. 1. These values are obtained from Equation IV-A-3. 1. 
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TABLE I \f-A~ 3. 





g d 


i/C_ 

lJ 


1/C Cj 

x> 


10 


0.5513 


* . 3 j 4 


o * 1 


12 


0.4823 


2.074 


6.33 


14 


0.4247 


2 . 3.4 


1 « C4 


15 


0.4026 

- 


2 . 484 


» V ■' 



Inspection of Figure IV-A-3.I shows the nterceuii on and V:_s 

limit cycle to occur at,, 

£0 ” 0# t’iv t ^ a f s e c 

and 

~ ^ • Jt 1 £ 5 r> 

The describing function method predicts a stable iia.-i . cyc^a > 
the following arguments. For an input signal grease r ». ha tn. the need 
zone v the relay characteristic is similar tc saturation and the erg 
ments given in Chapter I v Section I-A-3 can be applied to Figure 




IV-A-4 Analog Computer Simulation 



The analog computer simulation of the given nonlinear s/ctero 
shown in Figure IV-A-4.1. The coefficient pet settings vita iic.xo iced 
resistances and capacitances for real tir/ie and real magnitude scaling 

are given in Table IV-A-4.1. 



TABLE IV-A-4.1 



Pot 


Settings 


Associated Elements 


a l 


£ <, 

i 

100 


Fo r 
R , 


an i\psi 6cep cr E 
- 1 Xega k f1 - i 


^6 


<L 


= 1.0 


R 2 


= i ffeg, K :i •- 1 


FL. & 


a„ 

j 


a 1.0 


R, 

j 


= 0.1 Hcg s R_2 “ - 


Me* 


a 

4 


- 1.0 


R 

u 


-0.1 Keg 9 R =- 1 


K-.g 


a 5 


- 1.0 


fv 

J 


“ i fieg 5 ^ ~ 1 


^ f 


& 6 


= 0.5 


R e 


= J Meg 5 u f „ « x 


ui 


d 

/ 


= 1.0 


K 

/ 


- i Meg s C ... ~ a 

1 jL 


ut 


00 


= 1.0 


R. 

6 


- ] Keg, 4 f . - 1 


uf 


*9 


— 1,0 


S 


" - C f3 i 


ul 


*10 


= 1.0 


R 10 


=0.5 Mag ^ C f3 = 1 


. r 

Ofe A. 


a H 


- 0.5 


R n 


- 0.1 Keg s C - 1 

In- 


af 


a i2 


“ 0.4 


R i2 


~ ^ 8 i iVie Eu ^ f ~ ^ 


x *1 E 


3 


= 1.0 


R i3 


~ 0.1 C - ~ 1 

w ' O 


uf 


~14 


- 0.8 


R 14 


- C. : :;ei f c - 1 

r ^ 


uf 


a 15 


- 1.0 


R 15 


1 Heg * C f6 " * 


*Jt L. 


a rr 


C.05 (Initial 


For 


'Dead Zone of 3 Vol i 


£S 




Sett ing) 


Both Side of Zero 




1 

a 


= 0,05 


For 


p / o 1 1 s Re lay 0 u. t p . 


LiC 






Vol 


tage 




(xH 


itial Setting; 
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To accurately set the back-biasing pots ‘ ? ’a a si 

generator was connected to the input c£ the relay ton. 

An X-Y recorder was then calibrated to 2 volte per inch on each ex. r 

The Input and output of the relay simulator were then connected to the 

1 

respective X and Y axis terminals. The pots <v a n and * a ,J were tne n ad- 
justed until the X-Y recorder plot closely approximated the static 
characteristic of the chosen relay with dead zone. The static chaisw-tec - 
istic obtained by the X-Y recorder is shown in Figure T/-A-4.2* 

The system was first excited by a 10 volt step input signal and th 
system was found to build up into a limit cycle as shown in recorder trau 
IV-A-4.1. The system was then excited by a 20 volt step input signal and 
found to go into the same limit cycle as shown by recorder trace IV-A-4.2 
From these recorder traces the limit cycle was found to have an angular 
frequency of 0.455 rad/ sec and an amplitude of 14 . 3 volts. 

The phase portrait of the system is shown in Figure itf-A-4.3. 
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RECORDER MARK H 





Relay with Deac£ Zon<2 It Esrroit Cba'*, r ^ j 
with No Zeros a 20 Volt Step Xnp^t . 

SCALES • ■ 

P 

£ - 1 Volt/Lir.e £, - 0.3 Volt /i ' 

flirt - s Seecr-ds/'Cl- 



3 Volts/Li 



IV-A-5 



Comparison of the Four Methods its. Lint <_yui.es P- cd'.^n. ion for 
a Nonlinear System of Figure Iv-A. 



Method 


Frequency 
rad/ sec 


Arp 1 1 1 ude 
(/ones) 


Mitrovic 0 s 


0. 4607 


10.;': 


Root Locus 


0.4592 


1 3.01 


Describing Function 


0.455 


14.1 


Analog Simulation 


0.455 


14.3 
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IV-B Loop Transfer Function with One Zero. 

The block diagram of the system chosen is shown in Figure IV-B. 



% 



• 4 * 



T 




Figure IV-B 

Relay with Dead Zone in Error Channel Loop 
Transfer Function with One Zero 
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IV-B-1 Mitrovic® s Method. 



The characteristic equation of the given system is 



sV i2.rs 3 + 3$s a +(i& +rooA/)s + fooM - 0. 

(IV-B'1.1) 



Assuming the last two coefficients variable, the characteristic 

equation becomes, 

S 4 f i2,rs 2 +3E£ 2 +Bi S + Bo •= O (iv-b-i.2) 

where 



B = 500 N (IV-B-1.3) 

o 

and 

Bj = 16 + 500 N # (XV-B-1.4) 

From Appendix A, Mitrovic®s equations for B^ and are 

0) ] (IV- B- 1.5) 

and 



IS, =• 3%a)*<hC < 6)lr\7ifLO*<fa C&J+ (iv-b-i.6) 

Substitution of the values of the 0 functions from Appendix A for 
as 0 gives the parametric equations of the stability curve as, 

B o = 38 ujfl 2 + (IV-B-1.7) 

and 

B = 12.5 Ud^ 2 (IV-B-1.8) 

At the intersection of the M-point locus and the stability curve, 
equation IV* B- 1,3 equals Equation IV-B-1.7 and Equation IV-B-1.4 equals 
Equation XV“B~lo8, thus 

500 N = 38 60* 2 + 60^ 4 (IV-3-1.9) 

and 

16 + 500 N = 12.5 60 n 2 (IV-B-1.10) 
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Substituting Equation IV-B-1.9 into Equation IV-B-1.10 and combining 
terms yields 

r O . (iv-b-i.H) 

2 

Using the quadratic formula to solve Equation IV-B-1.11 for Ldr\ 
yields 

w A 2 - 

from which the frequency of the limit cycle is determined to be 

=5*11 rad/sec. 

2 

Substituting (j) n back into Equation IV-B-1.9 the "average 11 gain of 
the nonlinearity over a cycle is determined as 

N ~ 0.6208 

from which the approximate amplitude of the limit cycle is determined 
to be 

£ = 8.054 volts, 

max 

A graphical presentation of the intersection is sketched in Figure 
IV-B-1.1. 

Mitrovic®s Method predicts a stable limit cycle by applying the 
same arguments given in Section IV-A-1 to Figure IV-B-1.1. 



. 
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Graphical Solution for Relay with Dead Zone in Error 
Channel, Loop Transfer Function with One Zero 



I 
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IV-B-2 Root Locus Method 



Substituting S = 3(0 in the characteristic Equation IV-B-i.l 
and evaluating the higher powers of J gives, 

LQ /l ~MZ.ruj 3 -Z8'co :L +Jtl(>+5‘00N)u) + TO ON-O, (iv-b-2.1) 

Requiring that the real and imaginary parts of Equation IV-B-2.1 
go to zero independently yields the following equations, 

LO 4 -3 8 60 ^ -h^OOA/ =6 (IV-B-2.2) 

and 

^!U+S-OOH)lX> =0 (IV-B-2.3) 

Eliminating N between equations IV-B-2.2 and IV-B-2.3 gives, 

L0^-2rS(^ a -1C* (IV-B-2.4) 

2 

Using the quadratic formula to solve Equation IV-B-2.4 for 60 
gives, 

60 2 = 26.11 

from which the frequency of the limit cycle is determined to be 

U) = 5.11 rad/sec. 

2 

Substituting 60 back into Equation IV-B-2.2 gives, the "average" 
gain of the nonlinearity over a cycle as 

N = 0. 6208 

ave 

from which the approximate amplitude of the limit cycle is determined 
to be 

f - 8.054 volts, 

max 

A graphical presentation of the root locus solution is sketched 
in Figure IV-B-2.1. 

The root locus method predicts a stable limit cycle by applying the 
arguments given in Section IV-A-1 to Figure IV-B-2.1. 
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Root J,ocus Solution for Relay with Dead Zone in 



Error Channel Loop Transfer Function with One Zero 
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IV-B-3 Describing Function Method. 

The describing function for a relay with dead zone is given in 
Equation IV-A-3.1* The linear loop transfer function for the given 
system in frequency response form is 

(C( JlL l)- 3h 2SM^} l 

^ JUJ + 1 )(0 , 25 J CO f j)(0. 1 25\) u>+ 1.) 

Table IV-B-3* 1 gives the magnitude and phase of Equation IV-B-3. 1 for 
the various frequencies used to plot the linear loop transfer function on 
the gain-phase plane of Figure IV-B-3. 1. These values are obtained from 
A Bode diagram plot of Equation IV-B-3. 1. 



TABLE IV-B-3. 1 







/G?i cieg 


4.65 


6 


-176 


4.95 


5 


-178.5 


5.1 


4.6 


-180 


5.3 


4 


-181.5 


5.6 


3 


-184 


5.9 


2 


-187 



Table IV-B-3. 2 gives the values of 1/G^ which are plotted in Figure 
IV-B-3. 1. These values are determined from Equation IV-A-3.1. 

TABLE IV-B-3. 2 





s 


1/g d 


1/G d db 


8 


0.621 


1.61 


4.14 


8.9 


0.592 


1.69 


4.56 


9 


0.588 


1.7 


4.62 


10 


0.551 


1.814 


5.175 
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Inspection of Figure IV-B-3,1 shows the intersection and thus the 



limit cycle to occur at 



and 



CO =5,1 rad/sec. 



f = 8.95 Volts. 
v — max 

The describing function method predicts a stable limit cycle by 
applying the arguments given in Section IV-A-3 to Figure IV-B-3.1. 
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IV-B-4 Analog Computer Simulation. 

The analog computer simulation of the given system is shown in 
Figure IV-B-4. 1. The coefficient pot settings with associated resistances 
and capacitances for real time and real magnitude scaling are given in 
Table IV-B-4. 1. 



TABLE IV-B-4. 1. 



Pot Setting 


Associated Elements 


E i 


E 

For an input step of i volts 


a i = Too 


Rj = 1 Meg s R fl = 1 Meg 


a 2 = 1.0 


R 2 = 1 Meg, R fl = 1 Meg 


a^ - 1«0 


R 3 = 0.1 Meg, R 2 = 1 Meg 


a. ■ 1.0 
4 


R =0.1 Meg, R = 1 Meg 

4 f J 


a^ = 0.3 


R 5 = 0.1 Meg, C fl = 1 uf 


a, =0.4 
6 


R & = 0.1 Meg, C fl *= 1 uf 


a = 1.0 
7 


R^ = 0.1 Meg, = 1 Meg 


a 8 = 1.0 


Rg = 0.1 Meg, R^ = 1 Meg 


a 9 = 1.0 


R g =0.1 Meg, C 2 = 1 uf 


a io = °’ 5 


R 10 = 1 Meg, C f2 = 1 uf 


a u = 0.5 


R n = 0.5 Meg, C f3 = 1 uf 


a 12 = °' 8 


R 12 = 0.1 Meg, C f3 = 1 uf 


a 13 = l.o 


R 13 = 1 Meg, G = 1 u£ 


a = 0 « 0 5 (Initial 


For 5 volt dead zone either 


Setting) 


side of Zero 


a" = 0,05 (Initial 


For 5 Volt relay output Voltage 


Setting) 





To accurately set the back-biasing pots "a" and "a ", the same pro- 
cedure as given in Section IV-A-4 was followed. The static character- 
istic of the relay is the same as that shown in Figure IV-A-4. 2. 
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9-HDO 




The system was first excited with a 6 volt step input signal and 
was found to build up into the limit cycle shown in recorder trace 
IV-R-4 0 1 o The ssytem was then excited with a 12 volt step input signal 
and the system was found to go into the same limit cycle as shown in 
recorder trace XV~B-4 0 2 e From these two recorder traces the limit cycle 
was determined to have an angular frequency of 5,07 rad/sec and an 
amplitude of ShO volts 0 

The phase portrait of the system is shown in Figure XV-A-4*2. 
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IV-B-5 Comparison of the Results of the Four Methods for the System 



of Figure X^-B. 



Method 


Frequency 
rad/ seCo 


Amplitude 

(Volts) 


Mitro/ic 0 s 


5.11 


8.054 


Root Locus 


Soil 


8.054 


Describing Function 


Sol 


6.45 


Ana log Si itju I at ton 


5.07 


9.0 




CHAPTER V 



NONLINEARITY IN A FEEDBACK PATH 
V-A Saturation in Velocity Feedback Path. 

The system chosen to analyze is given in Figure V-A. 




V-A- 1 Mitrovic's Method. 

Assigning N as the instantaneous gain of the nonlinearity, the 
characteristic equation can be written as, 

S 4 + 7s 3 + 14s 2 + ( 8 + 120N)S + 30 = 0. (V-A-l.l) 

Assuming the last two coefficients variable, the characteristic 
equation is rewritten as, 

S 4 + 7s 3 + 14s 2 + B 1 S + B Q = 0 (V-A-1.2) 

where 

B t = 8 + 120N (V-A-1.3) 

and 

B q = 30 . (V-A- 1.4) 

From Appendix A, Mitrovic's equations for B q and B^ are, 
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(V-A-1.5) 



8 0 - -D4(A) n 2 <^ / ^) + 7a)n S ^iW 4UJ» 4 <0 3 (<ej7 

and 

G,= H &>„$,&) + 7 w,, a <£ 3 «J) + Oi^^cq) . (v-A-1.6) 

Substituting the values of the 0 functions from Appendix A for £ = 0 
gives the parametric equations of the stability curve as 

B q = 14 UJrl^ (V-A-1.7) 

and 

B x = 7L0H 2 ; (V-A-1.8) 

From equations V-A-1,3 and V-A-1.4 the M-point locus is seen to be 
a straight line segment parallel to the axis. At the intersections 
of the M~point locus with the stability curve, equation V-A-1.4 equals 
V-A-1.7 and equation V-A-1.3 equals equation V-A-1.8, Thus at the inter- 
sections 

4 ' L 

60^ 4 14 = 36 (V-A-1.9) 

and 

7 d 0 n ' 1 = 8 + 1X0 N * (v-A-1.10) 

Solving equation V-A-1,9 for gives 

L 0^= 2.6411 

and 

Lt)n = 11.3589, 

from which 

£0^ = 1.625 rad/sec 

and 

U)/\ 2 _ = 3.370 rad/sec. 



148 



Since the system is fourth order, the stability curve is a para- 

bala with its axis parallel to the B axis. Consequently, both values 

o 

of (jO ^ must be considered since the M-point locus, a straight line 
parallel to the B 1 axis, can intersect the parabola in two places. 

For = 1*625 rad/sec, the average gain of the nonlinearity 

over a cycle is found from equation V-A-1.10 to be 



N, = 0.0874, 
lave 

from which the approximate amplitude of signal X in Figure V-A is deter- 
mined to be 



lmax 



= 34.325 volts, 



For (A* 3 .370 rad/sec, the average gain of the nonlinearity over 
a cycle is determined to be 



N 2ave - 0.5958, 

from which the average amplitude of signal X is determined to be 

X 0 = 5.035 volts. 

2max 

However, the signal X is proportional to the output velocity and 
not the output angle. To obtain an approximate amplitude of the limit 
cycle in terms of output angle, 0^ , one must assume the signal X is 
relatively sinusoidal. With this assumption, the amplitude of the output 
angle can be obtained from 

X * 






&.£> to 



(V-A-l.ll) 



From equation V-A-l.ll, the approximate amplitudes of the output 
angle for both limit cycles are 

= 3.52 volts 

max 

and 




Qofcnax * °- 2M volts ' 
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The graphical presentation of the solution by Mitrovic's method 



is shown in Figure V-A-l.l. 




FIGURE Y-A-U 

GRAPHICAL Solution for saturation tn 
velocity fee ob A c k path 



Mitrovic’s method predicts that the limit cycle at 1.625 rad/sec 

is an unstable one and that the limit cycle at 00^ = 3.37 rad/sec is a 
stable limit cycle by the following arguments. For the stable limit cycle 
at = 3.37 the same arguments used in Chapter I, Section I-A-l can 

be applied to Figure V-A-l.l. For the unstable limit cycle assume Qji 
is such that the first peak of signal X causes the instantaneous M-point 
to be at M^. The system is then stable and the oscillations decrease. As 
the oscillations decrease N increases moving the M-point to the right. 
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Eventually, the "average" M-point will lie at the intersection of the 
stable limit cycle at £0^^= 3.37, Now assume that Qji is large enough 
to cause the first peak of signal X to be large enough to place the 
instantaneous M-point at . The system is unstable and the oscillations 
increase* As the oscillations increase, N decreases causing the M-point 
to move even farther to the left into the unstable region. Consequently 
the oscillations increase indefinitely. 
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V-A-2 Root Locus Method, 

Substituting S = J into equation V-A-l.l and evaluating the 
higher powers of ^ gives, 

Lb 4 -J7<x> 2 - J4uj z •H2 qn)0J+ 30 = 6, (v-a- 2. i) 

By requiring that the real and imaginary parts of equation V-A-2. 1 
go to zero independently, then 

CO 4 - 14 U i 2, +*30 = 0 (V-A-2. 2) 

and 

-07 60^ -f 3(9+ /2-OAJ ) to -o , (V-A-2. 3) 

Dividing equation V-A-2. 2 by JtO the above two equations are 

exactly the same as equations V-A-1.9 and V-A-1.10 developed in Mitrovic’s 

method. The solutions of these simultaneous equations are again given 

CO i ~ 1.625 rad/sec 

LQ 2 = 3.370 rad/sec 

N- = 0.0874 
lave 

N 0 = 0.5958 
2ave 



X. = 34.325 volts 
lmax 



and 



X 0 = 5.035 volts. 
2max 



Also applying equation V-A-l.ll gives 

©/> , =3.52 volts 

u lmax 



and 



S 0 2„,ax ’ °- 249 VOltS ' 

The graphical solution of the root locus solution is given in Figure 
V-A-2. 1. In the figure, the "open loop poles" are the poles of the tran- 
sfer function obtained when the position feedback loop is closed on the 
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) 



plant, first manipulating the 1.5 gain block out of the error channel. 




The root locus method predicts the limit cycle at 60j = 1.625 rad/sec 
will be an unstable limit cycle and the limit cycle at 3.37 rad/sec 

will be a stable one. The arguments of Chapter I, Section I-A-2 can be 
applied to the 3.37 rad/sec crossover to determine a stable limit, . 

cycle at that point. At the 60, = 1.625 rad/sec crossover, assume the 

input signal, , is large enough to cause the first peak of signal 

X to be large enough to place the instantaneous root at r^. The system 
is stable and the oscillations decrease. As the oscillations decrease, 

N increases causing the instantaneous closed loop root to move up the 
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locus toward the linear closed loop root. 

Eventually, the average root location over a cycle will lie at the 
stable limit cycle crossover. If is large enough to cause the first 

peak of signal X to be large enough to place the instantaneous root at 
r^ s then the system is unstable and the oscillations increase. As the 
oscillations increase, N decreases, driving the closed loop root farther 
into the right-half plane toward the open loop root. Consequently, the 
oscillations will increase indefinitely. 
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V-A-3 Describing Function Method. 

The describing function for saturation is given in Equation I-A-3.1 
of Chapter I. 

To obtain an equivalent loop transfer function, the gain of the 
error channel was first made unity by moving the gain block to the input 
channel and to the unity feedback path as shown in Figure V-A-3. 1. 

Reduction of the outer position feedback thus yields 



G (s) 



20 

' $ 4 + 7s 2 4 I4s 2 + ?s+ 30 • 



(V-A-3. 1) 




FIGURE Y-A-2./ 

A digital computer program was used to find the roots of the de- 
nominator of equation V-A-3. 1, which are, 
r = 0.2168 - jl. 3854 
r 2 = 0.2168 + j 1.3854 
r 3 = •- 3.7168 - j 1.2011 
r = - 3.7168 + j 1.2011. 

Thus the equivalent open loop transfer function becomes, 



G H- (s) — 



\ZOS 

(S-0.2I&8 ±J 1*3854 )(S F 3.I7C8 tJI.20il) 



(V-A-3. 2) 
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With S = Ju3 s the magnitude and angle of equation V“A-3 a 2 were 
calculated using a digital computer program with CO varying between 0*1 
rad/sec and 100 rad/sec „ The results of this program with angles in the 
neighborhood of “180° are given in Table V~A~3 0 l e These results are plot- 
ted on the gain-phase of Figure V-A-3„2 0 



TABLE V-A-3.1 



CD 


/fe'HOujJj J b 


/6'hOio) ^ 


1.60 


21.65 


-182.57 


1.61 


21.46 


-181.51 


1.62 


21.27 


-180.50 


1.63 


21.08 


-179.54 


1.64 


20.89 


-178.64 


1.65 


20.71 


-177.79 


3.30 


4.903 


-179.02 


3.35 


4.605 


-179.72 


3.36 


4.556 


-179.86 


3.37 


4.498 


-179.995 


3.38 


4.441 


-180.134 


3.39 


4.384 


-180.27 


3.40 


4.327 


-180.41 



Values of the X/G^ locus are given in the two areas of interest in 
Table V-A-3*2 U These values are calculated from Equation X-A-3.1* 
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TABLE V-A-3.2 



R 


Gp 


l/G o d & 


0.060 


0.07635 


22.34 


0.066 


0.08397 


21.52 


0.068 


0.08651 


21.26 


0.070 


0.08905 


21.01 


0.072 


0.09159 


20.53 


0.074 


0.09413 


20.53 


0.46 


0.5643 


4.969 


0.47 


0.5756 


4.798 


0.48 


0.5868 


4.630 


0.484 


0.5935. 


4.532 


0.488 


0.5957 


4.499 


0.490 


0.5979 


4.467 


0.50 


0.6090 


4.308 



By interpolation of Figure V-A-3.2 and the tables, two limit cycles 
are seen to exist at 

£<J| = 1.625 rad/sec 
/?, = 0.0686 

XjfrtAjf 43 - 7 volts 

and 



= 3.37 rad/sec 
= 0.488 

^2Wa.x = 6 * 15 volts * 
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To determine the approximate amplitude of the output angle* sinusoidal 
variations of the signal X must be assumed. Thus equation V-A-l.ll can be 
used to give* for the limit cycle at (j ) j = 1,625 rad/ sec 

lmax • 4 ' 48 volts 

and for the limit cycle at 3.37 rad/sec 

&0 2max = °' 304 toUs - 

The describing function method predicts a stable limit cycle at 
C0 A = 3 .37 rad/ sec and an unstable limit cycle at 6l)| = 1.625 rad/sec. 

At L0^= 3 ^ 37 the limit cycle is determined to be stable by applying the 
same reasoning given in Chapter I* Section I-A-3 to Figure V-A-3.2 using 
the signal X, At £0| “ 1,625 the limit cycle is determined to be unstable 
by the following reasoning. If the input* * is large enough to 

cause the first peak of signal X to be large R will be small and the in- 
stantaneous critical point will lie at say in Figure V-A-3.2, The 
system is unstable due to a negat ivephase margin and the oscillations in- 
crease. As the oscillations increase* R becomes even smaller driving the 
critical point farther up the -180° axis and increasing the negative phase 
margin. Thus the system is always unstable and the oscillations increase 
indefinitely. If is small enough such that the first peak of signal 

X is small enough* the instantaneous critical point will lie at say R^ of 
Figure V-A-3.2. The system is stable due to positive phase margin and the 
oscillations decrease. As the oscillations decrease R becomes smaller 
driving the critical point down the -180° axis and increasing the positive 
phase margin. Eventually the positive phase margin begins to decrease due to 
the shape of the Q O CO,) curve and the system reverts to the 

stable limit cycle at 60^= 3.37 rad/sec. 
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V-A-4 Analog Computor Simulation. 



The analog computor simulation of the given system is shown in 
Figure V-A-4. 1. The coefficient pot settings and the associated resist- 
ances and capacitances for real time and magnitude scaling are given in 
Table V-A-4. 1. 



TABLE V-A-4. 1 



Pot Setting 


Associated Elements 




For input of Ei Volts and Ke « 1.5 


0.75 Ei 




a l “ 100 


Pxj = 0.5 Meg, = 1 Meg 


a 2 = 0.75 


R 2 = 0.5 Meg, R fl = 1 Meg 


= 1.0 


R 3 = 1 Meg, R f2 = 1 Meg 


a. = 1.0 
4 


R^ = 1 Meg, R f2 = 1 Meg 


a $ = 1.0 


R 5 = 1 Meg, C fl = 1 uf 


a, = 1.0 
6 


Rg = 1 Meg, C fl = 1 uf 


a_ = 1.0 

i 


R ? = 0.5 Meg, C f2 = 1 uf 


a g = 1.0 


R g =0.5 Meg, C f2 = 1 uf 


a 9 = 1.0 


R g = 0.1 Meg, C f3 = 1 uf 


a 10 = °* 4 


R 10 = 0.1 Meg, C f3 = 1 uf 


a ll “ ! *° 


R j3 = 1 Meg, R f3 = 1 Meg 


a l2 = 


R l2 =0.1 Meg, R f ^ = 1 Meg 


a 13 “ 1 ’° 


R 13 = 1 Meg, R f5 = 1 Meg 


o 

IJ 


R 14 = 1 Me S> c f 4 = 1 uf 


a = 0.03 


For a 3 Volt Saturation 


Initial Setting 


Voltage 



To accurately set the back-bias pots "a", a sine wave generator 
was connected to the input of the nonlinearity simulation and the output 
monitored. Pots M a" were then adjusted until the output was limited to 
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Ke-U 




Analog Computer Simulation of Saturation in Velocity Feedback Path 



3 volts. The static characteristic curve is the same as that in Figure 
X~A.4 W 2, except the saturation voltage is 3 volts instead of 5 volts. 

A step input signal of 1 volt was first impressed on the system and 
the system was found to build up into the limit cycle shown in recorder 
trace V-A-4.1, A step input signal of 5.8 volts was then impressed and 
the system was found to go into the same limit cycle as shown in recorder 
trace V-A-4.2. As expanded portion of the stable limit cycle is shown in 
recorder trace V-A-4.3. 

A visual presentation of the unstable limit cycle was difficult to 
obtain. For a step input of 6.21 volts the system was only slightly 
damped before reverting to the stable limit cycle as shown in recorder 
trace V-A-4.4. For step input of 6.22 volts the system was slightly un- 
stable with the oscillations eventually increasing indefinitely as shown in 
recorder trace V-A-4.5. 

From the recorder traces, the stable limit cycle was determined to be, 

60 x = 3.35 rad sec 

X = 6.2 volts 

7 lmax 

, = 0.30 volts. 

Ob l ma x 

The unstable limit cycle was determined by averaging recorder traces 
V-A^4 tt 4 and V-A-4.5 with measurements taken on the first few cycles only, 

6U| = 1.618 rad/sec 

Xi^ x = 42 « 3 volts 

Qoi 4.29 volts. 

The phase portrait of the system is shown in Figure V-A-4.2. 



162 





A 2921 32 BRUSH IXSTBUMEMTS division of clevite corporation 




Saturation in Velocity Feedback Path, 1 Volt Step Input. 

* 

SCALES; £-0.1 Volts/Line £ - 0.1 Volts/Line X - 0.5 Volts/Line 

t . 

Y - 0.2 Volts/Line Time - 1 Second/Division 
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Recorder Trace V-A-4«2 

Saturation in Velocity Feedback Path, 5,8 Volt Step Input* 

# 

SCALES“£-0«5 Volts/Line £ -0,5 Volts/Line X - 5 Volts/Line 

I r 

Y - 0.2 Volts/Llne Time - 1 Second/Division 
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FI INSTRUMENTS division of cieviu corporation . cievei, 





Recorder Trace V-A~4*3 

Saturation in Velocity Feedback Path* Limit Cycle . 

SCALES: £ -0.02 Volts/Line £ - 0.1 Volts/Linfc ' X - 0.5 Volts/Line 
f - 0.2 Volts/Line Time - 0.2 Seconds/Division 
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Recorder Trace V-A-4.4 
Saturation in Velocity Feedback Path, 6.21 Volts Step Input. 

SCALES: £ - 0.5 Volts/Line £ - 0.5 Volts/Line X - 5 Volts/Line 

i - 0.2 Volts/Lme Time - <1 Second/Division 

— - - — — 1 6 6u_ 




Recorder Trace V“A~4 0 5 

Saturation in Velocity Feedback Path, 6.22 Volts Step input. 

SCALES t- £“-0.5 Volts /Line £ - 0.5 Volts/Line X - 5 Volts/Line 

i - 0.2 Volts/Line Time - 1 Second/Division 
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V-A-5 Comparison of the Four Methods in Predicting Limit Cycles of the 



System of Figure V~A. 



Method 


Angular 


Frequency 




Amp lit uce (7 c Its) 




r ad/ sec 




X max 


E max 




Stab 1 e 


. Unstable. . . 




.. ...Stable . 


- .Gast&blj* 


. Stable— 


ILasbahle. 


Mitrovic & 


3.37 


1.625 




5.03 


34.33 


0.249 


3.52 


Root Locus 


3,37 


1.625 




5.03 


34.33 


0.249 


3.52 


Describing 


3.37 


1.625 




6 , 15 


43.7 


0.304 


4.48 


Function 
















Analog 

Simulation 


3.35 


1.615 




6.2 


42.3 


0.30 


4.29 
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It 

V-B Saturation in Acceleration Feedback Path. 

The system chosen to analyze is shown in Figure V-B. 




Figure V-B 

Saturation in Acceleration Feedback Path 



V-B-l. Mitrovic's Method. 

Assigning N as the gain of the non-linearity the characteristic 
equation is: 



4 3 2 

s + 7s + (14 + 10N)s + 8s + 20 = 0 „ 

Assuming the s and s 2 coefficients to be variable the 
equation becomes 

4 3 2 

s + 7s- + B 2 s + B^ + 20 = 0 

where 



(V-B-l. 1) 
characteristic 

(V-B-l. 2) 

( 



and 




(V-B-l. 3) 



B 2 = 14 + 10 N. , 



(V-B-l. 4) 
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From Appendix A* Mitrovic°s equations for B. (1 and are 



a, 




2 Ofaltf) -7ui« 3 <f> l Cg)-u),fyl<i)] 



(V-B-1.5) 



and 



Bl-±z[ ^ / (-*£) t 7C0 k1 (V-B-1.6) 



Substitution of the 0 fur. sticks from Appendix A for {f =0 gives 



the parametric equations of the stability curve as 



and 



B i - * 1 COki 

_ 20 



(V-B-1.7) 



'1 " avS + (On 






(V-B-1.8) 

The shape of the stab! 1 icy curve on the versus E 2 plane Is not 
easily recognized. However, at the intersection of the stability curve 
and the M- point locus ^ equation v-B-1.7 must equal equation V-B-1.3 or 

7 ^ , >. 

from which the frequency of the Unit cycle is determined to be 

(X)f\ -* 1.06$ rad/sec. 

Also at the intersection equation V-B-1.8 must equal equation V-B-1.4 
or 

20 4 

14 - ION = ~-U -j- 00 n 1 0 

Od* 

Sub st itut icv of LO yields the average gain of the nonlinearity over 
a cycle of 

N ~ 0.^614 

ave 

from which the approximate amplitude of signal X is determined to be 



X = 10.82 volts 
- max 

To obtain an approximate amplitude of the angular oscillations, one 
must assume X to be sinusoidal. Then the amplitude of the error signal 
may be determined by 



or 




X fr\<x* 

os tof 



(V-B-1.9) 



S> = 18.94 volts, 
max 

The graphical solution by Mttrovic's method is sketched in Figure 

V-B-l.l. 




Acceleration feed back path 
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Mitrovica's method predicts an unstable limit cycle by the following 

reasoning* If the input signal is large enough such that the first 

peak of signal X is large enough, the instantaneous nonlinear gain will be 

small enough to place the M-point at M in Figure V-B-l.l the system is 

a 

unstable and the oscillations increase. As the oscillations increase N de- 
creases moving the M-point further into the unstable region. Thus the 
average M-point can never be in the stable region and the oscillations 
increase indefinitely. If the input signal is small enough such that the 
first peak of signal X is small enough, N will be relatively large placing 
the instantaneous M-point at M^ in figure V-B-l.l. The system is stable 
and the oscillations decrease. As the oscillations decrease N increases 
approaching one and moving the average M-point farther into the stable 
region. Thus the system is always stable and all oscillations eventually 
die out. 
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V-B-2„ Root Locus Method,, 



Substitution of S = JCO into the characteristic equation and 
evaluation of the higher powers of J yields 

<JU 4 -J 7Cu 3 -(I4 + ION)Go 2 + .} 8 60+ 2 0-0. (V-B-l.l) 

Requiring that both the real and imaginary parts of equation V-B-l.l 
go to zero independently 

cu 4 _ (if +ioa/)( J d 2 + 20 -O <v-b-i.2) 

and 

— v\7 60^ + 8 tu = 0 , (V-B-1.3) 

By dividing equation V-B-1.3 by JLD , the frequency of the limit cycle 
is obtained directly as 

60— lf^7 — l,0£>9 r<=id/ Seo . 

Substituting into equation V-B-1.2, the average gain of the non- 
linearity is obtained as 

N = 0.4614 
ave 

from which the approximate amplitude of signal X is obtained as 

X =10.82 volts, 

max 

To get the approximate amplitude of the error signal, X is assumed 

sinusoidal and equation V-B-1.9 is used, 

g = 18.94 volts, 

max 

The root locus solution is sketched in Figure V-B-2.1. Before 
sketching the figure, the unity feedback loop was first reduced and the 
poles of the reduced transfer function were used as the poles of the 
equivalent loop transfer function. 
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The root locus predicts an unstable limit cycle by applying the same . 
reasoning given in Section V-A-2 to Figure V-B-2.1. 
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V-B-3. Describing Function Method. 

The describing function for saturation is given in equation I-A-3.1 
of Chapter I. 

To obtain an equivalent loop transfer function the unity feedback 
loop was first reduced which gives, 




20 

S 4 -+7S 3 + 14 S 2 '4 8 S +2.0 



(V-B-3. 1) 



A digital cotnputor program was used to find the roots of the de- 
nominator of G 5 (s). 



r t = 0.0869 + j 1.201 
r 2 = 0.0869 - j 1.201 
r 3 = -3.5869 + j 0.9652 
r 4 = - 3.5689 - j 0.9652 

Thus the equivalent open loop transfer function becomes 
r'n,r\ - „ 'Q3 2 - 

m5J " (S-0.086<?±J |,20l)( S f 3. 



(V-B-3. 2) 



With 5 ~ J U> , the magnitude and angle of equation V-B-3. 2 were 
calculated by a digital computor program with iO varying between 0.1 
rad/sec and 100 rad/sec. The results of this program with angles in the 
neighborhood of -180° are given in Table V-B-3. 1. These results are also 
plotted on the gain-phase plane of Figure V-B-3. 1 as the linear loop trans- 
fer function. 
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TABLE V-B-3.1 



60 


1 6'Moal il= , 




1.05 


5.599 


-182.94 


1.06 


6.153 


-181.48 


1.07 


6.721 


-179.83 


1.08 


7.301 


-177.97 



The values of 1 /G d in the area of interest are given in Table V-B- 

3.2. 



TABLE V-B-3.2 



R 


g d l?3 


1/G n db 


0.30 


0.3762 


8.493 


0,35 


0.4364 


7.203 


0.37 


0.4601 


6.743 


0.372 


0.4625 


6.698 


0.374 


0.4648 


6.654 


0,38 


0.4719 


6.523 


- 0.40 


0.4954 


6.101 



By interpolation of Figure V-B-3.1 and the tables, the limit cycle 

is at 



60 = 1.609 rad/sec 

and 



from which 



R = 0.374 



X = 13.37 volts, 
max 
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To obtain an approximate amplitude of the error signal, X is 

assumed sinusoidal and equation V-B-1.9 used, 

£ = 23.4 volts, 

max 

The describing function method predicts an unstable limit cycle by 
applying the same reasoning given in Section V-A-3 to Figure V-B-3.1. 
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V-B-4. Analog Computer Simulation. 

A good simulation of the acceleration signal was difficult to obtain. 
The simulation that was finally decided upon was to use the velocity signal 
and approximately differentiate it by simulating s/ (s + 100). By using 
this transfer function an extra gain of 100 had to be simulated in the 
feedback path in order to maintain the D.C. gain constant of the path un- 
altered. Thus the simulated gain in the feedback path is 100 or 50 
instead of 0.5. 

The analog computor simulation is shown in Figure V-B-4. 1. The co- 
efficient pot settings and associated resistances and capacitances for 
real time and magnitude scaling are given in Table V-B-4. 1. 

The diode back-biasing pots were set in the same way described in 
Chapter I. The static characteristic curve is the same as Figure I-A-4.2. 

An input signal of 17 volts was first applied and the system was 
found to be stable but only slightly damped as shown in recorder trace 
V-B-4. 1. An input signal of 23 volts was then applied and the system was 
found to be unstable as shown in recorder trace V-B-4.2. The unstable 
limit cycle was found to occur for an input step of 19.8 volts and is 
shown in recorder trace V-B-4. 3. 

From recorder trace V-B-4. 3, the limit cycle is determined to have an 

amplitude of X = 10.6 volts and a frequency of £0 = 1.03 rad/sec. 
max 

The phase portrait of the system is shown in Figure V-B-4.2. 

The simulation of the acceleration signal is admittedly poor* especial- 
ly since the amplitude of X differs considerably from that predicted by 
the describing function. 
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Figure V-B-4.1. 

Analog computer Simulation of Saturation in Acceleration Feedback Path 



Table V-B-4.1 



Pot Setting 


Associated Elements 


Ei 

a i ~ 100 


For a step input of 
= 1 Meg, 


Ei Volts 
R fl = 1 Meg 


a 2 = 1.0 


R 2 = 1 Meg, 


R fl = 1 Meg 


a 3 = 1.0 


R 3 = 1 Meg, 


R f2 = 1 Meg 


a. = 1.0 
4 


R 4 = 1 Meg, 


R f2 = 1 Meg 


a 5 = 1.0 


R 5 = 1 Meg, 


C fl = 1 uf 


a 6 = 1.0 


R & = 1 Meg, 


C fi = 1 uf 


a 7 = 1.0 


= 0.5 Megj 


C f2 = 1 uf 


a 8 = 1 “° 


R g = 0.5 Meg, 


C f2 = 1 uf 


a g = 1.0 


R 9 = o.l Meg, 


C f3 = 1 uf 


a io = °" 4 


Rjq =0.1 Meg, 


C f3 = 1 uf 


a ll = 1,0 


R 13 =0.1 Meg, 


C.. = 0.1 uf 
f4 


a l2 = 1,0 


Rj 2 =0.1 Meg, 


C c/ = 0.1 uf 


= 1.0 


R 13 =0.1 Meg, 


R f3 = 1 Meg 


a l4 = K0 


R 14 = 0.1 Meg, 


R f3 = 1 Meg 


S 15 = °« 5 


R 15 = 0.1 Meg, 


R f4 = 1 Meg 


a 16 = 1,0 


R 16 = 1 Meg, 


R f5 = 1 Meg 


a i7 = 1,0 


R 17 = 1 Meg, 


C f5 = 1 uf 


a !8 = 1,0 


R lg = 1 Meg, 


R f6 = 1 Meg 


a - 0.05 Initial Setting 


For a saturation of 


5 volts 



Although the nonlinearity input X is not strictly sinusoidal, it is 
not too far from a sinusoid. The large discrepancy can not be accounted 
for solely by wave shape. However, if a root locus solution similar to 
that of Section V-B-2 is done on the system assuming 50s /(s + 100) in 
the inner feedback path the following results are obtained ; 
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60 = 1.0725 



When the results are compared with those of the root locus solution. 
Section V-B-2, it is seen that the extra pole in the feedback path causes 
a slight increase in frequency and a significant increase in the approxi- 
mate amplitude* Consequently, whatever it was in the actual analog simu- 
lation that caused the frequency to be less (1.03 rad/sec) than the pre- 
dicted value (1.069 rad/sec) also caused the amplitude to be considerably 
less than predicted by the describing function. 
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CLEVELAND. OHIO PRINTED IN U S A 




\ 

# 

SCALES 2 £ 2 volts/line £ - 2 volts/line 

y - 0.5 volts/line X - 1 volt/line 



Time - 1 second/division 




Recorder Trace V-B~4 C 1 

Saturation in Acceleration Feedback Path, 17 volt Step Input 
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RECORDER MARK H chart no. ra 2921 32 brush instruments 




SCALES: £ - 2 volts/line £ - 2 volts/line 

Y - 0.5 volts/line X - 1 volt/line 

Time - 1 second/division 







Recorder Trace V-B-4.2 

Saturation in Acceleration Feedback Path, 23 volt Step Input 



] 

» 
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9 \ 

SCALES 0 £ ~ 2 voits/line £ - 2 volts/line 

Y - 0o5 volts/line X - 1 volt/ine 

Time - 1 second/division 



r v/l 





Recorder Trace V-B-4,3 

| Limit Cycle* Saturation in Acceleration Feedback Path* 19,8 Volts Step Input 
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V-B-5* Comparison of the Four Methods in Predicting Limit Cycles of the 



System. of Figure V-B„ 



Method 


Frequency 
rad/ sec 


Amplitude (Volts) 
X-Max P, Max 


Mitrovie 9 s 


1.069 


10.82 


13.94 


Root Locus 


1.069 


10.82 


18.94 


Describing Function 


1.069 


13.37 


23.4 


Analog Simulation 


1.03 


10.6 


20.0 
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PART II 



A STUDY OF LIMIT CYCLES OF A SYSTEM WITH 
TWO GAIN VARIABLE NONLINEARITIES USING MITROVIC' S METHOD 

\ 

Introduction 

The author had no experience in the problem of a system with two 
non-linearities and consequently had no idea of what performance to expect. 
Thus a system was arbitrarily chosen and an analog computer study of the 
system was started. 

The system chosen is shown in Figure VI. This particular configura- 
tion was chosen primarily because each of the nonlinear gains appears alone 
in two of the coefficients of the characteristic equation, which is, 

s 4 + 7s 3 + 14s 2 + (8 + KK N 0 )s + KK N, =0 (Vl-a) 

t l el 




Figure VI 

System with Two Saturation Nonlinearities 
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Assuming the last two coefficients variable, the characteristic equa- 



tion becomes 

a fy 

§ + 7 a + 14s + B.s + B =0 (VX-b) 

1 o 

where 

B - KK N, (VX-C) 

o el 

and 

B 1 = 8 + KK fc N 2 (VX-d) 

The equations of the stability curve are the same as those in Chapter 
V„ Section V-A-l and are repeated below 

B q = 14 OJ^ 2 - U)„ 4 (Vl-e) 

and 

B 1 =7UJ n 2 (VI- f) 

This particular arrangement easily lends itself to presentation on the 
versus B^ plane in that 

lo The stability curve is independent of the gain K since the 

4 3 2 

stability curve only depends on the coefficients of the s s s , and s terms* 

2, Once K is chosen a second rectangular grid of constant and 

K can be drawn G 
e 

3* Once and are chosen then any instantaneous or average 
M~point is easily plotted using the grid* For example, if = 6 

and = 4 and by some means an M-point is determined by = 0*5 and = 
0*25 9 then this M-point is easily plotted at the intersection of the = 3 
line and the ~ 1 line of the secondary grid* 

In the study 9 three terms are used to describe certain parts of the 
system shown in Figure VI and need to be defined* These are: 

Plant “ that portion of the system represented by the linear trans- 
fer function ^ _ K 

} “ SCS+/XS + ZXS+4) 
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Linear System - the equivalent linear closed loop transfer 
function with both feedback loops reduced and no nonlinearities present, 
or the linear system described by the, characteristic equation vi-a with 

Nj and N 2 equal to 1.0. 

Linear M~point - that point on the B q versus B„ plane which is 
determined by equations VI-c and Vl-d with ^ and equal to 1.0. 
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Chapter VI 



ANALOG COMPUTER STUDY AND STATEMENT OF RESULTS 
Besides the poles of the plant, there are six parameters which directly 
or indirectly affect the two assumed variable coefficients of the character- 
istic equation* These are K^, ^sl* K 1 9 anc * ^ s 2 * As a starting point 

the plant gain, K, was made 20 and both saturation voltages were made equal 
to 3 volts* For these conditions the analog simulation is shown in Figure 
VI-1* The settings of the coefficient pots and their associated resistances 
and capacitances are given in Table VI- 1* 

The step input signal, or , was obtained by a pot reduction of a 
100 volt d« c a source through an isolating amplifier* To change various 
combinations of a^, a^, and were used, however, both channels had to 

be the same* To change K various combinations of a and R were used* 

t 1 j 1 j 

Because the diode resistance is not zero during conduction and pot 

loading effects of pots n a H and "a®", the back bias pots were readjusted 
slightly as the amplitude of the limit cycle changed* The readjustments 
were made by eye such that on the average the output signals of both satura- 
tion simulators were 3 volts* A typical output signal of a saturation simu- 
lator for a large input signal is sketched in Figure VI-2* 

First the stability curve was drawn on the B q versus plane and then 
a second grid of constant and lines was drawn as showa, in Figure VI-3* 
From this figure, it is apparent that one should get a good indication of the 
limit cycles possible by taking many linear starting points located by and 
K and observing the resulting limit cycles* 

The initial study took exactly this form, that is, a was chosen and 

several responses to a step input were observed for various values of K on 

t 



192 



TABLE VI- 1 



Pot Setting 

a = K for K 4, 1.0 

1 e e 

a = for K > 1.0 

1 10 e ' 

« K for K / 1.0 

2 e e ^ 

*2 - TO £ " K e > l -° 

a^ = 1.0 

a. = 1.0 
4 

a^ = 1.0 

a,. = 1.0 
6 

a « 1.0 

7 

a 8 * 1 *° 
a g “ 

*10 ■ 1,0 

a ll “ 0<> ^ 

a j2 - 1.0 

a X3 * K t f ° r K t 4 1.0 

a 13 - I fot K t > ‘-O 

8 14 ’ *-° 

8 15 ’ U0 

*16 * !.0 



Associated Elements 

R x = 1 Meg, R fl = 1 

Rj * 0.1 Meg s R fl = 1 

R 2 = 1 Meg s R fl = 1 

R 2 = 0.1 Meg s R fl = 1 

R 3 = 1 Megs R f2 - 1 

\ . 1 Meg, R f3 “ 1 

R 5 « 1 Megs R f3 = 1 

R 6 = 1 Meg, C fl = 1 

R ? - 1 Meg, C fl « 1 

R g =0.5 Meg, C f2 = 1 
R^ ~ 0.5 Megs C £2 - 1 
R^q = 0.1 Meg, C £3 = 1 
R u - 0.1 Meg, C f3 « 1 
R^2 38 1 Meg, ^f4 = ^ 

R 13 - 1 Meg R f5 = 1 

R 13 =0.1 Meg, R f5 = 1 

R 14 e 1 Meg, R f6 = 1 

R 15 = 1 Meg, C f4 = 1 

R 16 s 1 Meg, R f7 = 1 



Meg 

Meg 

Meg 

Meg 

Meg 

Meg 

Meg 

uf 

uf 

uf 

uf 

uf 

uf 

Meg 

Meg 

Meg 

Meg 

uf 

Meg 



a -- 0.03 Initial Setting 
a' = 0.03 Initial Setting 



For a 3 volt saturation level 
For a 3 volt saturation level 
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a horizontal traverse across the plane* In the study 9 was varied from 

« 0,5 to I s i0, K was varied from K s 0 o 5 to K » lQ„0 o About every 
e t t t 

other value of K 9 the input step*, E 9 was varied over several values to see 
t ' 1 

if the resulting limit cycle was dependent on the magnitude of the step in- 
put signal 0 The only limitation placed on E 0 was overloading of any one of 
the computer amplifiers; consequent ly & for certain linear starting points 
E was as large as 200 volts* The results of these computer runs are 

XiuSl K 

tabulated in Table I of Appendix Bo 

From this first portion of the analog computer study several things 
were evident * These were: 

1„ No limit cycles were obtained from a step input signal when 
the linear system M~point was inside the stable region on the versus 
plane no matter what the input signal magnitude* 

2* Limit cycles were only obtained for step input signals when the 
linear system M-point was outside the stable region and these limit cycles 
were independent of input signal magnitude and only dependent on the linear 
system M-point position*, 

3 0 The unstable region of the B versus Bo plane can be divided in- 

o I 

to three regions such that only one of the nonlinear elements saturate or both 
elements saturate for the linear system M-poxnt lying in a given region,. 
Several runs were made for various linear M-points inside the stable 
region and for various step input signal magnitudes 0 From the results of 
Section V-A of Chapter one might expect for linear M-points close to the 
left portion of the stability curve that the system could be driven unstable 
for large step input signals* However*, all linear M-points chosen showed 
some damping even for the linear point extremely close to the stability 
curve* A typical response for a linear M~point very close to the left 
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portion of the stability curve is shown in Recorder Trace VI -1 for a 
small step signal input and Recorder Trace VI-2 for a large step signal 
input, Typical responses for a linear M- point very close to the right 
side of the stability curve are shown in Recorder Trace VI-4 for a large 
step signal input. Typical responses for a linear M-point in the well 
damped portion of the stability region are shown in Recorder Trace VI-5 
for a small step signal input and Recorder Trace VI-6 for a large step signal 
input , 

Although it has been stated that no limit cycles were observed for a 
linear M-point in the stable region*, the system could be driven into a limit 
cycle for repeated disturbances in the input signal after the initial step 
input. This type of response was only noted for linear M-points with a very 
small damping ratio and close to the left portion near the maximum of the 
stability curve, A typical response of this type is shown in Recorder Trace 
VI-7 for extremely little damping. In this trace*, a limit cycle was achie- 
ved with only one disturbance after the initial step. Recorder Trace VI-3 
shows the same type of response but with more damping. In this trace three 
disturbances were required to achieve a limit cycle. 

To obtain these responses an initial step input signal was applied. 

Then with the system in a very oscillatory slightly damped condition the 
disturbance was applied. To apply the disturbance,; the input pet was given 
a small quick twist to approximate a step change in input signal. It was 
found that the disturbances had to be (1) generally small in amplitude such 
that the error channel did not appreciably saturate*, (2) applied at a defin- 
ite time during a cycle and in a specific direction. It was found that the 
maximum effect was achieved when a positive disturbance was applied when the 
error was passing through lero from minus to plus, Also*, the maximum effect 



was observed for a negative disturbance as the error was passing through 
zero from plus to minus * This type of system performance was only found 
late in the experimental phase of the study and was not investigated 
thoroughly. 

Recorder Traces VI-9 and VI- 10 show that the resulting limit cycle for 
a linear M-point in the unstable region is independent of the magnitude of 
the step input signal. These are typical responses * and at all points where 
the effect of the step input signal was investigated, the limit cycle was the 
same for different input magnitudes. Comparison of Recorder Trace VT-9 or 
VI-10 with Recorder Trace VI-11 illustrates the effect of the relative posi- 
tion of the linear M-point and the stability curve on the resulting limit 
cycle. Recorder Trace VI-11 shows a limit cycle of higher frequency and much 
lower amplitudes than the one in Trace VI-10, 

The most striking result of the initial computer study was that the 

unstable region of the B versus B- plane can be divided into three sub- 

o 1 

regions. These subregions are shown on the coefficient plane of Figure 
VI-3. For the linear M-point in Region I, only the error channel saturates 
and the limit cycle can be predicted from the intersection of the co= 
ordinate of the linear M-point and the stability curve. For the linear 
M-point in Region II both channels saturate, but the resulting limit cycle 
frequency always indicated an intersection of an average M-point locus with 
the stability curve at a point to the left of the maximum point of the 
stability curve. For the linear M-point in Region III, only the tacho- 
meter channel saturates and the limit cycle can be predicted from the inter- 
section of the coordinate of the linear M-point and stability curve. 

The dividing line which proceeds upward and to the left from the sta- 
bility curve will be referred to as dividing line No. 1, The dividing line 
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which proceeds almost horizontally to the right from the top of the 
stability curve will be referred to as dividing line No, 2. 

The approximate position of dividing line No, 1 was obtained from the 
results tabulated in Appendix B. To fix more accurately the position of 
the dividing line* more runs were taken on the analog computer simulation 
by choosing a and varying K^. A point on the dividing line was then 
determined when the amplitude of signal X 2 in the resulting limit cycle was 
equal to 3 volts* the saturation voltage. The results are tabulated in 
Table 2 of Appendix B. Typical changes in the limit cycle as was in- 
creased across No. 1 dividing line are shown in Recorder Traces VI-12 and 
VI-13, The main feature to note from these two traces is the slight in- 
crease in the amplitude of signal from just under 3 volts to just over 
3 volts. 

The approximate position of dividing line No, 2 was obtained from the 
results of Table I in Appendix B. To fix more accurately the position of 

the dividing line* more computer runs were made by fixing K and varying 

t 

The results are tabulated in Table 3 of Appendix B, Typical changes 

in the response as was increased above the dividing line are shown in 

Recorder Traces VI-14 and VI-15. Recorder Trace VI-14 is a couple of 

oscillations of the limit cycle resulting from a linear M-point at K - 8,0 

and K =2.38. However, sufficient time was allowed to ensure that the 
e 9 

system did not exhibit the same response as shown in Recorder Trace VI- 15, 

Recorder Trace VI- 15 is a series of three sections taken out of the response 

for a linear M-point at K = 8,0 and K =2.39 to show the salient features 
K t e 

of the entire response. The time intervals between the sections are in- 
dicated on the trace. Section one shows the build up to an apparently 
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stable limit cycle with only the tach channel saturating. However, closer 
inspection reveals that the amplitudes of signals X^ and are increasing. 
Section 2 starts at approximately the same amplitudes as occur at the end 
of Section 1, but quickly builds up to the final stable limit cycle as shown 
in Section 3 at an expanded time scale. The salient features of these two 
recorder traces are as follows: 

1. The sudden build up from an apparent limit cycle to a final 
stable limit cycle for a linear M-point just above dividing line No. 2. 

2. The large difference in signal magnitudes between the limit 
cycles of the two traces and the difference in their frequencies. 

3. The change in frequency of the oscillations during the build 
up to the final limit cycle in Recorder Trace VI-15. This point is not 
very noticeable on Trace VI-15 but it is noticeable on traces with an ex- 
panded time scale using dividers to measure the period. 

4„ The frequency and amplitude of Trace VI- 14 is very nearly the 
frequency and amplitude of the apparent limit cycle of Trace VI-15. 

Since the foregoing results allow one to predict the limit cycles for 
a large number of linear system M-points for 3 volt saturations, the next 
logical step was to see if the dividing lines are independent of saturation 
level. Thus both error channel and tach channel saturation voltages were 
adjusted to 10 volts. (pots "a" and "a 1 " equal to 0.1). The same analog 
computer procedures used to determine the dividing lines for - 3 

volts were again used. The resulting limit cycles for the linear starting 
points are tabulated in Tables 4 and 5 of Appendix B. A comparison of the 
dividing lines for the 3 volt and 10 volt saturation levels is given below 
I in Tables VI-2 and VI-3. 
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TABLE VI -2 



No, 1 Dividing Line 


10 volt Saturation Level 




3 volt Saturation Level 


K 


K 




K 


K 


e 


t 




e 


t 


2.1 


0.875 




2.1 


0.92 


2.2 


0.855 




-- 




-- 






2.3 


0.89 


2.4 


0.84 




-- 




2.5 


0.83 




2.5 


0.85 


3.0 


0.795 




3.0 


0.83 


4.0 


0.785 




4.0 


0.81 


5.0 


0.775 




5.0 


0.80 


6.0 


0.77 




6.0 


0.785 



TABLE VI- 3 



No, 2 Dividing Line 


10 volt Saturation Level 


3 Volt Saturation Level 


K 


K 


K 


K 


e 


t 


e 


t 


2.445 


2.5 




-- 


-- 


__ 


2.425 


3.0 


2.415 


4.0 


2.415 


3.0 


— 


-- 


2.405 


5.0 


2.395 


6.0 


2.395 


6.0 


2.385 


8.0 


2.385 


8.0 


2.335 


15.0 


2.335 


15.0 
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Tables VI-2 and VI-3 indicate that both dividing lines are essential- 
ly independent of saturation level, although dividing line No. 1 is more 
doubtful than dividing line No. 2. In order to definitely show this 
point, more saturation voltages will have to be used and more sophisticated 
measurement procedures should be used. 

Practically all of the data given in Appendix B is graphically dis- 
played in Figure VI-4 through VI-9. Figure VI-4 is a plot of the limit 
cycle amplitude of signal a s is varied with fixed. Figure VI-5 

is a plot of the limit cycle amplitude of signal for like variations. 
Figure VI-6 is a plot of the limit cycle amplitude of signal as is 
varied with fixed. Figure VI-7 is a plot of the amplitude of signal 
X^ for like variations of K^. Figure VI-8 is a plot of the limit cycle 
frequency as is varied with fixed. Figure VI-9 is a plot of the 
limit cycle frequency as is varied with fixed. 

These six figures essentially show the results as already stated. By 
correlation of these figures with the location of the linear M-point in 
Figure VI-3, one can see the regions in which only one element saturates. 

For example, the upper right portion of Figure VI-8 shows Region III of 
Figure VI- 3 in which one would expect the frequency to remain constant as 

is increased. Also, the transition across dividing line No. 2 is very ap- 
parent in all the figures because of the large discontinuity in the curve. 
Some curves for relatively small and abruptly stop and then start up 
again in a different portion of the graph. This is due to the variable 
being inside the region of stable linear M-points in Figure VI-3, thus no 
limit cycle. An important result to note which has not yet been pointed out 
is that for greater than 3 and K greater than 3, the frequency of the 
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resulting limit cycle is essentially independent of variations in either 

K or K . 
e t 

Figure VI-10 shows various instantaneous and "average" quantities plotted 
against the number of elapsed oscillation peaks for a linear M-point just 
above No. 2 dividing line. A better abscissa name would be n number of 
elapsed half cycles 1 * although both are related to the elapsed time by the 
characteristics of the system. The quantities plotted in the figure are 
defined below: 

B^min = The minimum instantaneous value during a half cycle 

of the s^ coefficient calculated by equation Vl-d 

using the minimum instantaneous gain . . N 0 . 

2min 2min 

occurs when the signal is maximum during the half 
cycle. 

The "average** value of the s^ coefficient over a half 

cycle as calculated from equation Vl-d using the "average" 

gain of the nonlinearity l^ave* The describing function 

was used to relate N 0 to the maximum value of X 0 

2ave 2 

during the half cycle. 

The minimum instantaneous value of the s° coefficient 
calculated by equation Vr-c using the minimum instan- 
taneous gain N, , . N- , occurs when the signal X- 
is maximum during the half cycle. 

The "average" value of the s° coefficient calculated 
from equation Vl-d using the "average" gain of the non- 
linearity N^ ave . The describing function was used to 

relate N, to the maximum value of X, during the half 
lave 1 



lave 



omin 



oave 



cycle. 
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0J = fc The average frequency of the oscillation over a half 
a v e 

cycle obtained by measuring the half-period of the 
oscillations. 

Justification for using the describing function in the manner indicated 
above will be indicated in Chapter IX. 

If the '’average" M-point is assumed to indicate the stability of the 
system over the half cycle, then it is interesting to note that the "aver- 
age" M-point, as defined by B„ and B , moves through a portion of 

lave oave 

the stable region on a constant B line. It is. also noted that under limit 

o 

cycle conditions at the extreme right side of Figure VI-10, the values of 

B^ and B q calculated using in Mltrovic f s stability curve equations 

are approximately equal to the end values of B- and B curves which 
rr J n lave oave 

were calculated using the describing function and the amplitude of the limit 

cycle. A comparison is made below. 

From Stability Curve From Figure VI- 10 

. - C. B- = 30.9 B- = 31.4 

1 lave 

B =* 43. 2 • B = 42.8 

o oave 




figure ye -2. 
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Recorder Trace VI- 1 

Linear M-Point Close to Left Side of Stability Curve* 

Ke = 2.4, = 1.7, 1,2 Volts Input Step 

SCALES : 

=0.2 Volts/line Y^ - 0.2 volts/line X 2 - 0.2 volts/line 

Y 2 - 0.2 volts/line Time - 5 Seconds/Division 
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Recorder Trace VI-2 

Linear M-Point Close to Left Side of Stability Curve, 



= 2.4 



K = 1.7, 62.5 Volts Input Step 
SCALES : 

- 10 Volts/Line - 0.2 Volts/Line * 0.5 Volts/Line 

-0.05 Volts/Line -0.05 Volts/Line -0.1 Volts/Line 

-0.05 Volts/Line 

Y~ - 0.2 Volts/Line 
-0.1 Volts/Line 
-0.05 Volts/Line 

Time - 5 Seconds/Division 
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Linear M-point Close to Right Side of Stability Curve 
= 1.0, = 3.9, 2.6 Volts Input Step. 

SCALES : 

X -0.2 Volts/Line Y 1 -0.2 Volts/Line Y ? -0.2 

1 -0.01 -0.01 Volts/Line 

X 2 -O. 2 Volts/Line 

Time - 5 Seconds/Division 



Volts/ 

Line 
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